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Objectives:  (List  the  objectives  of  the  research  effort.  This  may  be  omitted  if  there  has 
been  no  change.  State  new  or  revised  objectives  if  they  have  changed  and  why.)  No 
change  to  the  objectives. 


Status  of  effort:  (Brief,  under  200  words,  of  progress  towards  achieving  the  research 
objectives.) 

This  DURINT  project  produced  high  impact  technical  accomplishments.  The  notable 
accomplishments  include  building  (i)  new  nanostructures,  (ii)  study  of  optoelectronic  and 
magnetic  properties,  and  (iii)fabrication  of  polymer  based  non-linear  based  photonic 
crystals  and  photonic  circuitry,  and  (iv)  fabrication  of  photodetctor  as  well  as  diffraction 
grating  devices.  A  colloidal  route  was  devised  to  produce  nano-heterostructures  by 
integrating  plasmonic  (Au),  magnetic  (Fe304)  and  semiconducting  (PbS  or  PbSe) 
moieties  into  a  single  platform.  New  ferromagnetic  (low  temperature)  InMnP 
nanoparticles  were  synthesized  and  characterized  by  far  IR  and  Raman  spectroscopies 
For  photodetector  devices,  an  organic/inorganic  polymeric  nanocomposite  consisting  of 
poly-N-vinyl  carbazole(PVK),  lead  selenide  quantum  dots(QDs)  as  photosensitizer  and 
organic  semiconductor  pentacene  as  a  conductivity  booster,  is  constructed.  The  device 
exhibits  dramatic  enhancement  of  infrared  photocurrent.  Another  device  consisting  of 
PVK,  PbSe  QDs  conjugated  to  carbon  nanotube  also  shows  a  significantly  enhanced 
photodector  efficiency.  In  further  work  on  the  device  front,  we  demonstrated  through 
holographic  photopolymerization,  the  formation  of  3-D  photonic  crystal  of  exceptionally 
long  range  order  and  high  efficiency  1-D  reflection  gratings.  Additionally,  we  have 
employed  two-photon  polymerization  technique  for  defect  engineering  in  3-D  photonic 
crystals.  Moreover,  a  dramatic  increase  in  third  harmonic  efficiency  was  achieved. 

In  addition  to  a  large  number  of  publications,  the  project  also  produced  a  number  of 
inventions. 
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1.  Objectives 


This  project  has  the  goal  of  providing  the  military  with  a  hybrid  technology 
involving  integration  of  molecular  electronics,  nanoelectronics,  nanophonotics  and  nano¬ 
optoelectronics.  It  is  a  materials  development  program  to  enable  appropriate  technology 
development.  The  objective  of  the  proposed  research  is  a  systematic  design,  at  both 
molecular  and  nanoscopic  levels,  leading  to  generation-after-next  materials  and 
structures,  in  which  electron  or  photon-mediated  processes  can  be  judiciously  controlled, 
not  only  to  enhance  the  performance  of  a  particular  device  but  also  to  lead  to  new 
capabilities.  The  objectives  for  materials  development  involves:  (i)  developing  novel 
molecular  and  supramolecular  structures  that  will  allow  a  judicious  tailoring  of  electronic 
and  photonic  interactions  and  dynamics  leading  to  control  of  channels  for  electron  and 
photon  (excitation)  transfer,  and  (ii)  developing  new  self-assembly,  field  driven  assembly 
and  processing  techniques  to  enable  production  of  periodic,  aperiodic  and  engineering 
architectures  on  the  nanoscale,  which  can  lead  to  co-operative  amplification,  synergism, 
new  manifestations  and  broad  spectral  responses.  In  this  program,  we  have  emphasized 
the  materials  development  aspects  of  nanotechnology  by  pre -tailoring  the  optoelectronic 
properties  materials  for  futuristic  applications. 


2,  Approaches 

The  various  tasks  of  the  project  require  a  multidisciplinary  approach  involving 
modeling,  synthesis,  process  and  characterization.  The  necessary  comprehensive 
approach  involves  the  following: 

•  Rational  design  of  molecular  structure  and  precursors  based  on  structure-property 
relationships. 

•  Novel  chemical  synthetic  approach  using  solution  chemistry  to  produce 
molecular,  dendritic  and  inorganic  semiconductor  nanostructures. 

•  New  processing  schemes  to  produce  inorganic/organic  hybrid  nanocomposites 
with  direct  interfaces. 

•  Processing  of  polymer  based  self-assembled  ordered  nanostructures  and  photonic 
crystals  for  nonlinear  optical  functions. 

•  Use  of  bioinspired  novel  supramolecular  structures  as  templates  for  periodic  or 
aperiodic  self-assembling  of  electronic,  photonic  and  magnetic  units. 

•  Characterization  of  nanostructures  by  TEM,  SEM,  AFM  and  confocal  microscopy 
as  well  as  by  various  linear  and  nonlinear  optical  spectroscopies. 

•  Characterization  of  various  electronic,  photonic  and  magnetic  functions. 

•  Interactive  feedback  between  design,  synthesis  processing  and  characterization. 

Our  approach  also  places  a  very  strong  emphasis  on  utilizing  the  expertise  of  various 
team  members  in  a  complementary  way.  A  strong  collaboration  with  the  DoD  laboratory 
and  interactive  feedback  with  them  ensures  that  the  goals  of  the  project  are  fully 
compatible  with  the  needs  of  the  military. 
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Accomplishments: 

Significant  accomplishments  were  made  under  this  DURINT  program  in  the  tasks 
originally  proposed.  A  major  emphasis  was  placed  on  working  as  a  highly  collaborative 
team  and  interacting  with  the  Air  Force  Laboratory.  Many  major  accomplishments  of 
this  project  would  not  have  been  realized  without  this  collaborative  effort.  Constant 
communication,  exchange  of  information,  and  interactive  feedback  were  maintained 
through  a  number  of  means  including  exchange  visits  by  researchers,  focused  review 
workshops  and  organization  of  symposia  at  professional  societies. 

A  large  number  of  publications  and  invited  lectures  by  the  DURINT  team 
exemplify  the  great  deal  of  research  accomplished  under  this  program  and  its  broad 
impact  on  science  and  technology,  internationally.  In  addition,  a  number  of  patents 
resulting  from  these  accomplishments  have  been  filed. 

I.  Two-Photon  Energy  Harvesting  Dendritic  Systems 

Collaboration:  P.  N.  Prasad,  University  at  Buffalo; 

J.M.J.  Frechet,  UC  Berkeley; 

L.-S.  Tan,  AFRL,  Dayton 

Two-photon  absorption  excited  intramolecular  energy  transfer  and  light 
harvesting  effects  are  demonstrated  for  the  first  time  in  two  types  of  Dendritic  Systems. 
These  systems  contain  both  an  antenna  structure  that  strongly  absorb  two-photon  energy 
at  -800  nm  in  the  near-IR  region,  and  transfer  energy  efficiently  to  a  core  moiety  where 
the  energy  can  be  more  effectively  utilized  for  performing  significantly  enhanced 
photonic  function  such  as  lasing  or  photochemistry.  These  systems  are  bio-inspired  from 
nature’s  photosynthetic  systems  utilizing  chlorophyll  antenna,  except  that  our  system 
utilizing  nonlinear  two-photon  absorption  that  permits  photons  in  the  near-IR  to  IR  to  be 
absorbed  and  transforms  efficiently  into  excitation  normally  created  by  a  higher  energy 
photon  in  the  visible.  In  our  dendritic  system  the  core  moiety  emits  efficiently  at  —590 
nm.  The  two-photon  energy  absorbed  by  the  antenna  structure  is  resonantly  transferred 
to  the  core,  where  the  emission  intensity  of  the  latter  is  enhanced  by  8,  20,  and  34  times 
for  the  three  dendritic  systems,  respectively.  This  type  of  novel  dendrimer  structures  and 
intramolecular  energy  transfer  mechanisms  can  be  useful  for  multi-photon  pumped 
frequency-upconversion  lasing,  as  well  as  for  highly-sensitive  optical  sensors  working  at 
IR  spectral  ranges. 

Another  dendritic  system  harvesting  two-photon  excitation  for  photoactivation  at 
the  core  involves  a  two-photon  chromophore  synthesized  at  AFRL,  Dayton  (Tan)  and  a 
catalytic  phorphyrin  core  moiety.  The  IR  photons  are  strongly  absorbed  by  the  two- 
photon  absorbing  chromophores,  and  the  energy  is  efficiently  transferred  to  the  porphyrin 
moiety  which  generates  highly  reactive  singlet  oxygen.  This  reactive  singlet  oxygen  is 
very  important  for  a  number  of  useful  photo  oxidation  processes,  including  the  biological 
applications  in  photodynamic  therapy. 
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II.  Band-Gap  Engineering  and  characterization  of  inorganic  quantum  dots 

Collaboration:  P.  N.  Prasad,  A.  Cartwright,  B.  Weinstein,  B.  McCombe  and  H.  Luo, 

University  at  Buffalo 

New  Precursor  Routes  for  Quantum  Dots: 

Prasad’s  group  has  been  involved  in  the  development  of  novel  and  efficient  organic  based 
precursor  routes  for  band  gap  engineering  of  quantum  dots  and  core-shell  structures, 
whereby  quantum  dots  absorbing  in  different  spectral  range  and  dispersible  in  common 
organic  solvents  can  be  produced. 

There  has  been  a  well  focused  work  on  the  synthesis  of  highly  luminescent 
quantum  dots  using  a  new  precursor  route  involving  a  fast  reaction  of  In(Carboxylate)3 
with  P(SiMe3)3  in  octadecene.  These  reactions  are  reproducible  and  result  in  relatively 
monodospersed  crops  (12-13  %  relative  standard  deviation)  of  quantum  dots  in  as  little  as 
2  hours.  This  procedure  does  not  need  the  addition  of  an  external  surfactant,  because  the 
carboxylate  group  generated  in  situ  bind  on  the  surface  of  the  nanoparticles.  It  should  be 
noted  that  in  the  conventional  route  using  TOPO  as  the  surfactant,  InP  quantum  dots  took 
up  to  6  days  to  prepare,  after  which  a  multi  step  work-up  was  required  to  obtain  relatively 
monodisperse  crops. 

Along  with  the  preparation  of  InP  quantum  dots,  we  have  developed  a  new 
method  for  the  preparation  of  CdSe  quantum  dots  and  CdS  quantum  dots.  This  method 
results  in  high  quality  particles  that  are  highly  monodisperse  (less  then  8%  relative 
standard  deviation).  Included  in  this  publication  we  have  also  developed  a  one-pot 
synthesis  of  CdSe/CdS  quantum  dots  Core/Shell  quantum  dots  which  is  completed  in  less 
than  3  hours. 

A  considerable  amount  of  time  has  been  spent  on  the  synthesis  of  InP/II-VI 
core/shell  quantum  dots.  We  have  synthesized  InP/CdS,  InP/CdSe,  InP/ZnS  and 
InP/ZnSe  core-shell  quantum  dots.  All  of  these  QDs  have  been  successfully  synthesized 
and  characterized  by  IR,  UV-vis,  photoluminescence  and  photoluminescence  emission 
spectroscopy. 

Highly  efficient  III-V  quantum  dots  for  Imaging: 

Prasad’s  group  has  focused  on  the  development  faster  synthetic  routes  of  highly 
luminescent  quantum  dots  and  their  applications.  This  novel  method  that  uses  no 
externally  added  surfactant  was  used  for  the  synthesis  of  highly  monodispersed 
hydrophillic  InP-ZnS  nanocrystals  and  were  used  as  luminescence  probes  for  imaging. 
Hydrophobic  InP-ZnS  nanocrystals  are  prepared  by  the  new  method  that  yields  high 
quality,  luminescent  core-shell  nanocrystals  within  6-8  hours  of  total  reaction  time.  Then 
by  carefully  manipulating  the  surface  of  these  passivated  nanocrystals,  aqueous 
dispersions  of  folate  conjugated  nanocrystals  (folate-QDs)  with  high  photostability  are 
prepared.  As  an  illustration  of  imaging,  we  demonstrated,  using  confocal  microscopy, 
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the  receptor  mediated  delivery  of  folic  acid  conjugated  quantum  dots  into  folate  receptor 
positive  cell  lines  such  as  KB  cells.  These  folate-QDs  tend  to  accumulate  in 
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Figure  1.  Confocal  images  showing  fluorescence  of  QD-FA  in  KB  cells  and  a  localized  PL 
emission  spectrum.  Ffere  green  channel  shows  the  transmission  images,  while  the  intensity 
coded  (red  to  white )  channel  shows  the  fluorescence. 


multivescicular  bodies  of  KB  cells  after  6hours  of  incubation  (Fig.  1).  Receptor  mediated 
delivery  was  confirmed  by  comparing  with  the  uptake  of  these  particles  in  folate  receptor 
negative  cell  lines  such  as  A549.  Efficient  two-photon  excitation  of  these  quantum  dot 
particles  was  also  demonstrated.  The  use  of  these  InP-ZnS  nanoparticles  and  their 


Wavelength  (nm) 


Figure  2.  Absorption  and  emission  spectra  of  InP(QD)  and  InP-ZnS(QDs)  dispersed  in 
chloroform.  The  inset  in  emission  spectra,  shows  the  surface  state  emission  from  bare 
InP  particles,  which  is  absent  in  case  of  ZnS  coated  core-shell  particles. 
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efficient  two-photon  excitation  can  be  potentially  useful  for  deep  tissue  imaging  for 
future  in-vivo  studies 

Phonon  Spectroscopy  of  Quantum  Dots 

A  collaborative  effort  with  Professors  Bruce  McCombe  and  Bernard  Weinstein  has 
produced  significant  progress  in  understanding  the  vibrational  modes  of  the  InP/II-VI 
core-shell  nanoparticles.  This  understanding  was  achieved  through  a  combination  of 
variable  temperature  IR  and  Raman  measurements  and  detailed  theoretical  modeling, 
making  use  of  a  recently  published  dielectric  continuum  model.  These  results  allow  us  to 
determine  the  ratio  of  the  radii  of  the  core  to  the  shell,  and  also  in  some  cases  to 
determine  that  the  chemical  processing  did  not  yield  just  core-shell  nanoparticles,  but 
rather  a  mixture  of  InP  nanoparticles,  II-VI  nanoparticles  and  core-shell  nanoparticles.  A 
particularly  interesting  result  is  an  apparent  resonant  Raman  effect  that  tunes  dramatically 
with  temperature  in  the  InP/ZnSe  core  shell  structures.  (Fig.  3).  At  room  temperature  a 
very  sharp  line  associated  with  a  ZnSe  surface  mode  (cob)  is  seen;  below  100  K  this  mode 
disappears  and  is  replaced  by  two  other  sharp  modes,  one  associated  with  the  InP  core 
(coc)  and  one  with  another  surface  mode  in  the  ZnSe  (©a).  This  dramatic  behavior, 
illustrating  tunable  vibronic  couplings  to  spatially  different  vibrational  modes,  is  a  result 
of  the  interesting  band  structure  of  the  core  shell  structure  with  a  small  valence  band 
offset.  These  results  indicate  that  the  vibronic  coupling  can  be  tailored  in  appropriate 
core  shell  structures,  and  this  in  turn  means  that  there  may  be  interesting  photonic 
applications. 


Wavenumber  (cm'1) 

Fig  3.  Raman  spectra  of  InP  quantum  dots 
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Phonons  play  an  important  role  in  nanophotonics,  as  they  control  the  dynamics  of 
photon  induced  processes.  We  continued  analysis  of  the  vibrational  modes  of  the  InP/II- 
VI  core-shell  nanoparticles  determined  by  IR  and  Raman  studies,  and  initiated  studies  of 
II-VI  core  shell  nanoparticles.  These  results  have  allowed  us  to  determine  the  ratio  of  the 
radii  of  the  core  to  the  shell,  and  also  in  some  cases  to  determine  that  the  chemical 
processing  did  not  yield  just  core-shell  nanoparticles,  but  rather  a  mixture  of  InP 
nanoparticles,  II-VI  nanoparticles  and  core-shell  nanoparticles. 

We  have  also  studied  the  vibrational  modes  of  II-VI/II-VI  core-shell  nanoparticles 
by  IR  spectroscopy  and  find  good  agreement  with  previous  studies.  We  have 
investigated  the  dependence  of  the  absorption  strength  of  these  well-characterized  core¬ 
shell  structures  on  the  filling  fraction  in  the  Csl  matrix.  Results  show  systematic 
increases  in  absorption  with  filling  fraction;  for  a  filling  fraction  of  about  6%,  the  In-P 
vibrational  bands  are  overabsorbed,  demonstrating  the  possibility  of  using  IR 
spectroscopy  to  investigate  local  Mn-P  modes  of  the  interesting  ferromagnetic 
nanoparticle  system  (InMn)P  synthesized  by  Y.  Sahoo  of  Prasad’s  group.  The  Mn- 
related  modes  are  expected  to  be  much  weaker 
than  the  In-P  modes  due  to  the  low  concentration 
(typically  a  few  percent)  of  Mn  replacing  In;  for 
filling  fractions  of  6%  or  higher,  these  modes 
should  nevertheless  be  observable.  This  could 
serve  as  an  important  characterization  tool  in 
understanding  the  magnetic  behavior  in  this 
system. 

Finally,  we  have  initiated  magnet- 
photoluminescence  studies  of  InSb  (a  narrow  gap 
semiconductor)  nanoparticles  produced  by  D. 

Lucey  of  Prasad’s  group.  These  nanoparticles 
luminesce  in  the  IR  and  are  much  more  amenable 
to  tuning  of  the  emission  wavelength  by  particle 
size  due  to  the  small  fundamental  energy  gap.  An 
example  of  PL  measurements  at  low  temperature  is  shown  in  Fig.  4.  The  main  peak 
appears  to  be  due  to  the  22  nm  nanoparticles,  and  the  low-energy  peak  is  associated  with 
large  clusters  of  nanoparticles  giving  rise  to  quasi-bulk  emission. 

III.  Plasmonic  Nanoparticles 

Collaboration:  P.  N.  Prasad,  M.T.  Swihart 
University  at  Buffalo 
Younan  Xia,  University  of  Washington 

Metallic  nanoparticles  and  nanostructures  form  a  major  thrust  area  of  nanophotonics, 
called  plasmonics.  They  produce  highly  confined  and  significantly  enhanced  optical 
responses.  We  thus  placed  a  significant  effort  on  making  multifunctional  plasmonic 
structures 
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Figure  4:  PL  spectrum  at  4.2  K  of  InSb 
nanoparticles  embedded  in  a  Csl  matrix.  Data 
were  taken  with  a  BOMEM  DA-3  FTIR 
spectrometer. 
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Plasmonics  Nanoparticles 


Here  our  focus  is  to  prepared  metallic  nanostructures  which  can  be  used  to  enhance 
photonic  response.  Gold  nanoparticles  have  been  prepared  using  bifunctional  ligands  of 
the  type  X-R-SH  (X=  -COOH,  -OH,  -NH2)  so  that  the  surface  is  terminated  with  these 
functionalities.  The  ligands  bind  to  the  gold  surface  through  their  thiol  (-SH  end).  The 
coverage  and  binding  strength  of  the  surfactant  on  the  gold  surface  have  been 
investigated  using  thermogravimetric  analysis  and  FTIR  spectroscopy.  The  nanoparticles 
appear  to  have  excess  surfactants  on  their  surface,  but  these  do  not  necessarily  form  a 
bilayer  structure.  The  surfactant  conformation  allows  enough  room  for  the  particles  to  be 
bonded  to  other  surfaces,  such  as  those  of  another  nanoparticle,  even  with  some  excess 
surfactant  present  on  their  surface. 

We  synthesized  core-shell  structures  consisting  of  monodisperse  polystyrene  latex 
nanospheres  as  cores  and  gold  nanoparticles  as  shells.  Use  of  polystyrene  spheres  as  the 
core  in  these  structures  is  advantageous  because  they  are  readily  available  commercially 
in  a  wide  range  of  sizes,  and  with  dyes  or  other  molecules  doped  into  them.  Gold 
nanoparticles,  ranging  in  size  from  1  to  20  nm,  are  prepared  by  reduction  of  a  gold 
precursor  with  reducing  agents  such  as  sodium  citrate,  sodium  borohydride  or 
tetrakishydroxymethylphosphonium  chloride  (THPC).  Carboxylate  terminated 
polystyrene  spheres  are  functionalized  with  2-aminoethanethiol  hydrochloride  (AET), 
which  forms  a  peptide  bond  with  carboxylic  acid  groups  on  their  surface,  resulting  in  a 
thiol-terminated  surface.  Gold  nanoparticles  then  bind  to  the  thiol  groups  to  provide  up  to 
about  50%  coverage  of  the  surface.  These  nanoparticles  serve  as  seeds  for  growth  of  a 
continuous  gold  shell  by  reduction  of  additional  gold  precursor.  The  shell  thickness  and 
the  roughness  can  be  controlled  by  the  size  of  the  nanoparticle  seeds  as  well  as  by  the 
process  of  their  growth  into  a  continuous  shell.  By  variation  of  the  relative  sizes  of  the 
latex  core  and  the  thickness  of  the  gold  overlayer,  the  plasmon  resonance  of  the  nanoshell 
can  be  tuned  to  specific  wavelengths  across  the  visible  and  infrared  range  of  the 
electromagnetic  spectrum,  for  applications  ranging  from  the  construction  of  photonic 
crystals  to  biophotonics.  The  position  and  width  of  the  plasmon  resonance  extinction 
peak  are  well  predicted  by  extended  Mie  scattering  theory. 


Nanoshells  of  gold  and  silver  on  Dielectric  cores: 

We  present  a  systematic  study  of  the  plasmon  resonance  absorbtion  spectra  of  gold 
nanoshells  and  silver  nanoshells  on  the  surfaces  of  monodispersed  polystyrene 
microspheres  of  different  sizes  as  well  as  of  silver  nanoshells  on  free-standing  gold 
nanoparticles.  Gold  nanoshells  were  grown  on  polystyrene  microspheres  with  diameters 
of  188  and  543  nm.  The  commercially  available,  initially  carboxylate-terminated 
polystyrene  spheres  were  reacted  with  2-aminoethanethiol  hydrochloride  (AET)  to  yield 
thiol-terminated  microspheres  whereon  preformed  gold  nanoparticles  were  anchored  and 
grown  further  to  render  a  complete  shell.  Silver  nanoshells  were  grown  on  polystyrene 
cores  with  diameters  of  188,  296,  and  543  nm.  The  growth  of  silver  shell  required 
seeding  the  polystyrene  surface  first  with  gold  nanoparticles.  Reduction  of  silver  nitrate 
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onto  these  gold-decorated  microspheres  resulted  in  increasing  coverage  of  silver  on  the 
polystyrene  core.  The  nanoshells  were  characterized  using  transmission  electron 
microscopy  (TEM),  scanning  electron  microscopy  (SEM)  and  UV-vis  spectroscopy.  By 
varying  the  core  size  of  the  polystyrene  particles  and  the  amount  of  silver  reduced  onto 
them,  the  surface  plasmon  resonance  of  the  nanoshell  could  be  tuned  across  the  visible 
and  the  near-infrared  regions  of  the  electromagnetic  spectrum.  Gold  core  -  silver  shell 
nanoparticles  were  synthesized  by  reducing  silver  nitrate  onto  free  citrate-gold 
nanoparticles.  These  core-shell  particles  formed  chain-like  structures  with  a  plasmon 
resonance  absorption  that  could  also  be  tuned  across  the  visible  spectrum. 

Gold  nanorods: 

Seeded  growth  of  gold  nanorods(NRs)  has  been  accomplished  in  a  micellar  medium. 
Cetyl  trimethylamoniumbromide  (CTAB)  forms  micelles  where  the  growth  of  rod 
shaped  gold  nanoparticles  occurs.  AgN03  is  introduced  into  the  growth  solution  to 
enhance  the  population  of  the  NRs.  The  roles  of  non  ionic  surfactants  such  as  Tween  and 
Triton  have  been  examined  and  ascertained  that  with  their  increasing  perscentages  in  the 
growth  solution,  the  length  of  the  NRs  increases  monotonically  to  a  critical  limit,  after 
which  it  decreases  again.  The  growth  pattern  of  the  NRs  fits  into  the  soft  template  model 
rendered  by  the  mixture  of  CTAB  and  non  ionic  surfactants. 

Gold  Seeded  Semiconductor  -Magnetite  heterostructures: 

PbSe  nanostructures  are  potentially  useful  for  telecommunication  and  sensing 
applications  in  near  IR.  While  nanoparticles  and  nanowires  of  IV-VI  semiconductors 
such  as  PbSe  have  been  prepared  in  the  past,  nanorods  of  these  materials  with  controlled 
length  and  diameter  have  not  previously  been  reported.  Here,  PbSe  nanorods  with 
controlled  lengths  and  diameters  were  synthesized  using  Au  and  Au-Fe304  black-eyed 
pea-like  hybrid  nanoparticles  as  catalytic  seeds.  PbSe  nanostructures  with  different 
morphologies  (dots,  rods  and  stars)  have  been  prepared  by  changing  the  precursor  to  seed 
ratio. 

IV.  Silicon  Nanoparticles 

Collaboration:  P.  N.  Prasad,  M.  T.  Swihart 
University  at  Buffalo 

Nanostructured  silicon  dispersed  in  plastics  offers  another  opportunity  for  hybrid 
electroluminescent  devices.  Optically  clear  and  stable  dispersions  of  brightly 
photoluminescent  Si  nanoparticles  were  prepared  under  the  direction  of  Professor 
Swihart,  by  covalent  attachment  of  alkenoic  compounds  to  the  particles.  Quenching  of 
photoluminescence  by  ethylamine,  diethylamine,  triethylamine,  pyrazine,  and  piperazine 
was  investigated.  The  photoluminescence  was  quenched  by  the  action  of  these 
nitrogeneous  species,  but  in  some  cases  could  be  partially  restored  by  the  addition  of 
trifluoroacetic  acid.  The  extent  of  restoration  of  photoluminescence,  after  equilibrium  is 
reached,  was  independent  of  the  sequence  of  addition  of  the  amine  and  the  acid.  The 
photoluminescence  quenching  and  recovery  are  influenced  by  a  combination  of  basicity, 
polarity,  and  steric  factors  of  the  quencher  molecules.  The  quenching  and  subsequent 
restoration  occurs  gradually  at  room  temperature  and  it  takes  several  minutes  to  reach 
equilibrium. 
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V.  Nanocomposites  for  efficient  optoelectronics  and  photorefractivity 

Collaboration:  A.N.  Cartwright,  P.N.  Prasad,  University  at  Buffalo; 

M.  Durstock,  AFRL,  Dayton 


Hybrid  organiciinorganic  materials  have  emerged  as  a  novel  class  of  electronic 
and  optoelectronic  media  for  a  number  of  potential  technological  applications,  including 
solar  energy  conversion,  photodetectors,  electroluminescent  and  photorefractives.  They 
offer  broadband  tunability  of  the  wavelength  response,  increased  photogeneration 
efficiency  and  improved  charge  collection  efficiency.  We  have  already  shown  that 
semiconductor  nanoparticles  are  efficient  photosenisitizers  for  reduced  electric  field 
photogeneration  in  hybrid  materials.  Now,  we  have  also  shown  that  photoconductive 
composites  obtained  by  doping  of  a  host  polymer  with  semiconductor  nanoparticles 
exhibit  an  increase  in  the  effective  mobility  of  charge  carriers.  This  will  be  of  importance 
for  more  effective  operation  of  various  devices  including  broadband  solar  cells.  We 
studied  the  mobility  of  holes  in  a  model  nanoparticle-sensitized  hybrid  organiciinorganic 
system  consisting  of  poly(N-vinylcarbazole)  (PVK)  doped  with  quantum  dots  of 
cadmium  sulfide.  We  found  typical  features  of  the  charge  carrier  transport  like  strong 
dependence  of  the  carrier  mobility  on  the  field  and  temperature  as  well  as  on  sample 
thickness  which  is  a  signature  of  dispersive  transport.  However,  significant  enhancement 
of  the  effective  carrier  mobility  is  also  noticed  with  the  increase  of  nanoparticle 
concentration,  although  it  was  still  well  below  the  percolation  limit.  A  preliminary 
explanation  of  this  phenomenon  assumes  an  active  role  played  by  the  nanoparticles  in  the 
transport.  A  simple  theoretical  model  based  on  time-  and  mean-free-path  dependent 
mobility  was  proposed  to  account  for  this  surprising  result.  The  magnitude  of  the 
enhancement  is  found  to  be  consistent  with  this  model. 

In  order  to  produce  polymer:quantum-dot  nano-composites  providing  a  broad 
spectral  coverage  to  harvest  photons  over  entire  spectral  range,  we  have  focused  on  two 
areas  of  synthesis  and  processing: 

•  Development  of  novel  and  efficient  organic  based  precursor  routes  for  band  gap 
engineering  of  quantum  dots  and  core-shell  structures,  whereby  quantum  dots 
absorbing  in  different  spectral  range  and  dispersible  in  common  organic  solvents 
can  be  produced.  We  have  made  GaP,  InP,  silicon,  as  well  as  core-shell  structures 
InP/ZnS,  InP/CdS,  InP/CdSe  of  different  sizes,  producing  optical  absorption  and 
emission  over  a  broad  spectral  range. 

•  Processing  as  nanocomposites  in  an  efficiently  hole  transporting  polymer  to  create 
a  direct  interface  between  the  quantum  dots  and  the  polymer  for  facilitating 
charge  transfer  and  subsequent  migration.  For  this  purpose  we  are  experimenting 
with  different  capping  organic  ligands  which  can  thermally  be  driven-off  from  the 
surface  of  the  quantum  dots.  Another  approach  is  to  use  a  conjugated  structure  as 
the  capping  agent  whereby  ji-electrons  facilitate  the  charge-transfer  and 
subsequent  migration. 

Two  provisional  patents  have  been  filed  on  the  novel  synthetic  approach  to  band 
gap  engineering. 
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Photoconductivity  and  Photorefractivity  at  1.3  pm  Using  Quantum  Dot:  Polymer 
Nanocomposites: 


Due  to  the  relatively  large  Bohr  radius  (20  nm)  of  the  exciton,  PbS  quantum  dots  offer 
unique  access  to  the  regime  of  extreme  quantum  confinement  and  this  property  makes 
them  a  suitable  choice  to  photosensitize  photoreffactive  (PR)  polymeric  composites  at  the 
technologically  important  wavelengths  of  1.31  and  1.55  pm.  PbS  nanocrystals  were 
produced  by  reacting  the  precursors  via  a  hot  colloidal  synthetic  route.  High  quality, 
uniform  sized  nanocrystals  (-6-8  nm)  with  highly  crystalline  phases  were  obtained  (Fig 

5).  _ 
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Fig  5.  TEM  image  of  the  PbS  nanocrytals 


Fig  6.Enhancecf  pfioiosensitivity  5ph  of  the 
nanocrystalline  PbS  compared  to  previously 
reported  samples 


The  photorefractive  samples  of  composition  PVK:ECZ:DEANST:PbS  (48:15:35:2  wt  %) 
showed  greatly  enhanced  photoconductivity  (almost  an  order  of  magnitude)  at  the 
operational  wavelength  of  1 .34  pm  compared  to  previous  results  (Fig  3).  Pronounced  PR 
two-beam  coupling  effect  at  the  operation  wavelength  was  observed,  leading  to  very  high 
optical  gains  (Fig.  4).  A  response  time  faster  than  other  hybrid  photorefractive 
nanocomposites  is  also  observed. 
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Fig  7.  Electric  field  dependence  of  the  TBC 
gain  coefficient  T  in  the 
PVK:ECZ:DEANST:PbS  nanocomposite  at 
1.34  pm.  The  highest  gain  recorded  was  184.5 
cm-1  at  92.6  V/pm. 


Photoconductivity  and  Photorefractivity  at  1.55  pm  Using  Quantum  Dot:  Polymer 
Nanocomposites: 

a.  Excellent  photoconductivity  and  narrow  luminescence  in  solution-processed  hybrid 
organic: inorganic  nanocomposites  at  optical  communication  wavelengths 
The  photo-response  of  a  polymer  composite  was  extended  to  the  1.55  pm 
telecommunications  band.  Size  tunable  absorption  of  the  semiconductor  quantum  dots 
enabled  us  to  successfully  sensitize  the  hybrid  devices  and  tune  their  spectral  response 
across  the  infrared  range  between  800  nm  and  2  pm.  The  nanocomposites  are 
photoactive  at  infrared  wavelengths  with  narrow  emission  bands,  tunable  with  the 
quantum  dot  sizes.  The  photoluminescence  spans  a  wide  spectral  range  in  the  infrared 
(Figure  8).  Efficient  harvesting  of  infrared  photo-generated  carriers  led  to  high 
photocurrents.  In  particular,  we  observe  for  the  first  time,  high  photoconductive  internal 
quantum  efficiency  (among  the  highest  reported  so  far  in  this  spectral  range)  in  these 
hybrid  devices  at  the  technologically  important  wavelength  of  1.55  pm  (Figure  9), 
achievable  by  low-power  continuous-wave  illumination,  without  the  necessity  of  lock-in 
techniques. 
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Figure  8.  Normalized  room-temperature 
photoluminescence  spectra  of  a  series  of 
polymer-NC  composites,  demonstrating  size- 
tunable  emission  in  the  IR. 


Figure  9.  Main  Panel:  photoconductive  quantum 
efficiency  as  a  function  of  device  bias  at  excitation 
wavelengths  of  1.34  pm  (open  circles)  and  1.55  pm 
(filled  circles).  Inset:  the  ratio  of  photo-  to  dark 
conductivity  for  devices  as  a  function  of  applied  voltage, 
at  excitation  wavelengths  of  1.34  pm  (circles)  and  1.55 
pm  (triangles). 


b.  Photorefractive  response  of  a  polymer  device  to  1.55um  optical  communication 

wavelength:  an  organic:inorganic  nanocomposite 

Using  colloidally-fabricated  PbSe  quantum  dots  photosensitization  of  a  polymeric 
nanocomposite  at  the  optical  communication  wavelength  of  1 .55pm  is  achieved.  Efficient 
generation  of  infrared  photons  lead  to  high  photoconductivity  and  a  significant  figure  of 
net  gain  co-efficient  (Figure  10).  A  steady  state  diffraction  efficiency  of  -40%  was 
obtained  in  the  dynamic  refractive-index  gratings  (Figure  11).  The  dynamics  of  the 
grating  build-up  was  extracted  by  fitting  a  biexponential  function  to  the  temporal 
evolution  of  diffraction  efficiency  that  is  correlated  to  the  growth  of  the  internal  space- 
charge  field.  The  modulation  amplitude  of  the  refractive  index  is  also  computed  from  the 
diffraction  efficiency  data.  The  first  observation  of  net  optical  gain  and  significant 
diffraction  efficiency  at  this  important  wavelength,  accomplished  with  low-power 
continuous-wave  laser  beams,  makes  these  hybrid  nanocrystal-sensitized  devices  a 
potential  choice  for  infrared  imaging  and  optical  communication  applications. 
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FigurelO.  Electric  field  dependence  of  the  TBC 
gain  coefficient,  r,  in  the  hybrid  nanocomposite 
at  1.55  pm.  The  highest  gain  recorded  was  86.4 
cm-1  at  59.4  V/an'1. 


Figure  11.  Main  panel:  Steady  state  diffraction 
efficiency,  rj,  as  a  function  of  the  external  electric 
field  in  the  nanocomposite  at  1.55pm  .  The  solid 
line  is  a  theoretical  fit  to  the  experimental  data  . 
Inset:  Transient  growth  of  the  DFWM  diffracted 
signal  at  an  applied  field  of  40  V/an  '. 


c.  Solution  processed  hybrid  polymeric  nanocomposite  for  infrared  photodetection 


This  task  was  taken  to  enhance  the 
photosensitivity.  An  organic/inorganic  polymeric 
nanocomposite,  incorporating  a  soluble  precursor  to 
the  semiconductor  pentacene,  poly-N-vinyl  carbazole 
and  lead  selenide  quantum  dots  as  photosensitizer,  was 
fabricated.  Thin-film  devices  were  spin  coated  from 
the  composite  solution.  Following  thermal  conversion 
of  the  precursor  to  pentacene,  the  devices  exhibit 
dramatic  enhancement  of  infrared  photocurrent(Figure 

12).  Efficient  generation  of  carriers  coupled  with  high 
photoconductance  results  in  the  highest  quantum 
efficiency  reported  to  date. 


App?i©<f  (V) 


Figure  12.  Photocurrent  density  as  a  function 
of  applied  voltage  in  devices  with  the  same 
proportion  of  PVK:  pentacene  (3:1)  but 
various  amounts  of  PbSe  nanocrystals  as 
indicated  in  the  legend. 


VI.  Electroluminescense  in  Quantum  Dot:  Polymer  Nanocomposites 

Collaboration:  A.N.  Cartwright,  P.N.  Prasad,  University  at  Buffalo 

Our  DURINT  team  has  successfully  demonstrated  electroluminenscnce  in  PVK  doped 
with  quantum  dots  of  CdSe/CdS,  prepared  by  our  new  precursor  route.  A  three-layer 
device  structure  Al/PVK:quantum  dot/ITO  was  used.  The  results  on  electroluminescence 
(EL)  and  its  comparison  with  photoluminescence  (PL)  spectral  distribution  are  shown  in 
Figure  13.  This  redshift  from  PL  to  EL  is  commonly  observed  in  hybrid  devices. 
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Fig  13.  I-V  curve  (left)  for  electroluminescence  and  comparison  of  electroluminescence  and 
photoluminescence  spectra  (right)  for  the  CdSe/CdS  quantum  dot:PVK  polymer  nanocomposite. 


VII.  Polymer  Based  Magnetic  Nanocomposites  for  Electronics 


Collaboration:  P.  N.  Prasad,  H.Luo  University  at  Buffalo 
L.  Kempel,  Michigan  State  University 
M.  Alexander,  AFRL,  Dayton 

A  major  focus  has  also  been  to  utilize  magnetic  nanoparticles  containing 
polymers  for  radar  applications  in  electronics.  This  is  in  collaboration  with  the  Air  Force 
Research  Laboratory  at  Dayton.  Thus  we  have  utilized  nanochemistry  to  produce 
magnetic  nanoparticles  of  Fe304 ,  Ni,  and  Co.  We  have  also  tried  to  align  these  Fe304 
nanoparticles  in  the  size  range  of  8-12  nm  by  self-assembly  on  GaAs  substrates  in  the 
presence  of  strong  magnetic  fields.  A  long  range  ordering  is  observed  in  which  the 
particles  self-assemble  into  a  distribution  of  elongated  clusters  with  a  predominant 
orientation  lengthwise  along  the  field  direction.  Hysteresis  loops  measured  parallel  and 
perpendicular  to  the  alignment  direction  show  substantial  directional  dependence.  The 
coercive  fields  in  the  direction  parallel  to  the  alignment  field  are  larger  than  those 
perpendicular  to  it  by  57%  and  136%  at  100K  and  5K,  respectively.  A  broad  peak  is 
observed  in  magnetization  profiles  obtained  with  zero-  field-cooling. 

We  prepared  a  biocompatible  ferrofluid  containing  dye-functionalized  magnetite 
nanoparticles  that  can  serve  as  fluorescent  markers.  This  method  entails  the  surface 
functionalization  of  magnetite  nanoparticles  using  citric  acid  to  produce  a  stable  aqueous 
dispersion,  and  the  subsequent  binding  of  fluorescent  dyes  to  the  surface  of  the  particles. 
Several  ferrofluid  samples  were  prepared  and  characterized  using  Fourier  transform 
infrared  spectroscopy  (FTIR),  x-ray  photoelectron  spectroscopy  (XPS), 
thermogravimetric  analysis  (TGA),  BET  surface  area  analysis,  transmission  electron 
microscopy  (TEM),  and  SQUID  magnetometry.  In  addition,  confocal  fluorescence 
microscopy  was  used  to  study  the  response  of  the  fluorescent  nanoparticles  to  an  applied 
magnetic  field  and  their  uptake  by  cells  in  vitro.  Results  are  presented  on  the  distribution 
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of  particle  sizes,  the  fluorescent  and  magnetic  properties  of  the  nanoparticles,  and  the 
nature  of  their  surface  bonds.  Biocompatible  ferrofluids  with  fluorescent  nanoparticles 
enable  optical  tracking  of  basic  processes  at  the  cellular  level  combined  with 
magnetophoretic  manipulation,  and  should  be  of  substantial  value  to  researchers  engaged 
in  both  fundamental  and  applied  biomedical  research. 

For  magnetic  force  microscopy,  we  are  experimenting  with  the  system  itself,  so 
that  we  can  push  it  to  its  limit  to  study  single  particles.  During  this  process,  we  have 
tested  the  system  with  magnetic  films,  and  obtained 
several  important  findings.  For  example,  we  have  for 
the  first  time  discovered  interaction  between  well 
defined  magnetic  domains  and  their  subsequent 
movement.  Fig  14  shows  a  sequence  of  images  as  a 
function  of  temperature,  in  which  two  magnetic 
domains  attract  each  other,  merging  into  one, 
eventually  disappearing. 

We  have  been  able  to  prepare  stable  dispersions 
of  metallic  and  magnetic  nanoparticles  in  compatible 
solvents  so  as  to  blend  them  homogeneously  in 
polymeric  matrices.  For  example,  Au,  Pd,  Ni,  Fe304 
etc  nanoparticles  have  been  successfully  incorporated  in 
polymeric  matrices.  Frequency  dependences  of 
capacitive  impedance  and  relative  permittivity  of  these 
embedded  nanoparticles  have  already  been  measured. 

A  general  enhancement  of  permittivity  and  decrease  in 
impedance  are  detected. 


Fig  14.  MFM  images  of  two 
domains  (the  dark  lines), 
merging  together  as  the  sample 
is  cooled. 


VIII.  Photonic  Crystals  with  Enhanced  Nonlinear  Optical  Response 

Collaboration:  P.  N.  Prasad  and  A.  N.  Cartwright,  University  at  Buffalo; 

Y.  Xia,  University  of  Washington, 

E.  L.  Thomas,  M.I.T. 

Non-DURINT  Team  Collaboration:  R.  W.  Boyd,  University  of  Rochester;  J.  Haus, 

University  of  Dayton. 

Photonic  crystals  are  ordered  arrays  of  dielectric  nanostructures  which  have 
received  considerable  attention  worldwide.  Through  a  multidisciplinary  team  effort,  we 
are  focussing  on  novel  processing  methods  to  produce  nonlinear  photonic  crystals  where 
the  field  enhancement  as  well  as  anomalous  refractive  index  dispersion  near  the  band  gap 
can  be  used  to  produce  considerably  increased  nonlinear  response. 

The  team  effort  has  been  directed  towards  using  (i)  self-assembling  of  polymeric 
nanoparticles  on  a  template  (with  Y.  Xia)  and  holographic  lithography  (with  Prasad  and 
Thomas)  to  produce  a  three-dimensional  photonic  crystal  and  (ii)  a  block-copolymer 
approach  (with  E.  L.  Thomas)  to  produce  a  one-dimensional  photonic  crystal.  Another 
approach  utilizes  laser-directed  ordering  of  nanoparticles  to  produce  a  one-dimensional 
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photonic  crystal  (with  R.  Vaia  and  T.  Bunning  at  AFRL).  This  approach  and  holographic 
lithography  are  described  separately  in  the  section  of  “Photon ically  Aligned  Nanoarrays” 
and  “Multibeam  Holographic  Lithography,”  respectively. 

To  our  knowledge,  our  team’s  work  on  the  third-order  nonlinear  optical  processes 
is  the  first  reported  results  and  show  significant  enhancement  of  these  processes  in 
photonic  crystals.  Our  work  showed  that  amplification  in  the  two-photon  emission,  a 
third-order  nonlinear  optical  process,  appeared  at  the  edge  of  the  stop  gap  of  the  photonic 
crystal  which  was  produced  by  dye  solutions  infiltrating  a  polystyrene  photonic  crystal. 
The  amplification  was  interpreted  as  due  to  the  local  field  enhancements.  Our  results 
show  promises  that  a  photonic  crystal  medium  can  be  utilized  for  low  threshold  and 
spectrally  narrowed  two-photon  lasing. 

Third-harmonic  generation  (THG)  is  a  very  useful  technique  that  can  convert  the 
coherent  output  of  infrared  lasers  to  shorter  wavelengths  in  the  visible  and  near 
ultraviolet.  However,  the  smallness  of  the  third-order  optical  susceptibility  and  the  strong 
natural  dispersion  of  the  refractive  index  of  most  materials  have  prevented  the  practical 
utilization  of  one-step  third-harmonic  generation.  Therefore,  in  practice,  one  employs  a 
cascaded  two-step  process  to  produce  the  third  harmonic  with  high  conversion  efficiency. 

Having  perfect  phase  matching  of  the  pump  and  generated  signal  is  essential  for 
efficient  nonlinear  frequency  conversion.  One  way  to  achieve  this  condition  for  third- 
harmonic  generation  is  to  use  an  anomalous  dispersion  region  where  the  refractive  index 
of  the  medium  decreases  with  the  optical  frequency.  In  bulk  media,  anomalous  dispersion 
is  typically  accompanied  by  extremely  high  absorption  which  prevents  the  useful 
implementation  of  this  idea.  However,  in  a  photonic  crystal  anomalous  dispersion 
created  by  the  periodic  structure  is  not  accompanied  by  loss  and  strong  third-harmonic 
generation  can  be  obtained. 

We  observed  a  dramatic  enhancement  in  the  efficiency  of  third-harmonic 
generation  using  a  three-dimensional  (3-D)  polystyrene  photonic  crystal  medium.  When 
an  ultrashort  laser  pulse  with  the  proper  wavelength  passes  through  a  3-D  polystyrene-air 
photonic  crystal,  an  enhancement  (~25X)  of  single-step  third-harmonic  generation  is 
observed  at  the  short-wavelength  edge  of  the  photonic  bandgap.  The  origin  of  the 
enhancement  is  due  to  the  phase  matching  provided  by  the  periodic  structure  of  the 
photonic  crystals. 

We  have  also  utilized  a  one-dimensional  photonic  crystal  fabricated  by  laser- 
directed  ordering  of  liquid  crystal  nanodroplets,  to  produce  dramatic  enhancement  of 
third-harmonic  generation  near  the  stop-gap  at  the  fundamental  wavelengths.  Our 
theoretical  model  predicts  that  the  third-harmonic  enhancement  was  due  to  both  the  phase 
matching  and  the  pump-field  localization  within  the  structure. 


Moreover,  we  have  demonstrated  electrically  switchable  third-harmonic  generation  in 
low  refractive  index  contrast  one-dimensional  photonic  crystals  pumped  by  a  near- 
infrared  laser  beam.  The  observed  phenomenon  of  switching  is  described  separately  in 
the  section  of  “Electrically  Switchable  Nonlinear  Photonic  Crystals.” 
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Even  though  polymers,  used  in  our  experiments,  possess  only  a  weak  third-order 
nonlinear  susceptibility,  we  observed  a  strong  third-harmonic  generated  beam.  These 
circumstances  suggest  that  a  photonic  crystal  composed  of  a  conjugated  polymer,  which 
possesses  a  much  larger  third-order  nonlinearity,  could  produce  more  efficient  frequency 
conversion,  because  strong  nonlinear  material  properties  and  phase-matching  conditions 
can  be  provided  separately.  Even  thin  layers  of  such  photonic  crystal  medium  could  be 
used  to  fabricate  infrared  visualizing  cards  that  operate  at  a  practically  important 
wavelength  of  1.5  /urn .  Once  the  technology  of  fabrication  is  developed  to  produce 
photonic  crystals  of  high  optical  quality  and  of  appreciable  thickness  (a  few  millimeters), 
highly  efficient  frequency  conversion  of  the  output  of  coherent  infrared  light  sources  to 
the  short-wave  length  range  of  visible  light  and  the  ultra  violet  will  be  possible. 


Electrically  Switchable  Nonlinear  Photonic  Crystals 

One-dimensional  photonic  crystals,  which  are  formed  by  light-driven  photonic  alignment 
of  liquid  crystal  droplets  in  a  polymer  host,  can  provide  a  medium  for  electrically 
switchable  third-harmonic  generation.  The  strong  enhancement  was  observed  in  the 
reflection  as  well  as  in  the  transmission  geometries;  in  both  cases  at  the  high-frequency 
band  edge.  The  observed  phenomenon  of  enhancement  was  explained  theoretically  with  a 
coupled-mode  model.  The  enhancement  of  the  third-harmonic  generation  disappeared 
when  an  electric  field  was  applied,  revealing  the  electrical  switchability  of  the  third- 
harmonic  generation  (Fig  15).  The  strong  switching  properties  were  supported  in  the 
crystal  up  to  a  modulation  frequency  of  10  kHz  but  started  to  deteriorate  at  frequencies 
higher  than  3kHz. 


Fig.  15.  Intensity  of  the  third-harmonic  signal  measured  in  transmission  geometry  and  plotted  as  a 
function  of  time  (applied  voltage).  The  lower  curve  indicates  how  the  applied  voltage  was 
changing. 
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Photonic  Crystal  Defect  Engineering 

The  existence  of  periodic  variation  in  the  dielectric  constant  in  one,  two  or  three 
dimensions  may  cause  the  formation  of  a  photonic  bandgap.  This  bandgap  forbids  the 
propagation  of  light  in  certain  range  of  frequencies,  through  the  photonic  crystal  and 
gives  rise  to  a  range  of  interesting  effects  not  achievable  in  homogeneous  media.  By 
addition  or  removal  of  some  dielectric  material  in  the  photonic  crystal  structure  one  can 
introduce  defects  inside  the  crystal  which  disrupt  the  discrete  symmetry  of  the  lattice  and 
create  localized  modes  within  the  photonic  bandgap.  As  a  result,  when  light  of  frequency 
within  the  bandgap  propagates  into  the  defect,  it  becomes  confined  within  the  defect,  as 
the  crystal  does  not  allow  extended  states  at  this  frequency. 

Recently  we  have  conducted  systematic  studies  of  the  two-photon  polymerization 
technique  for  defect  engineering  in  3D  polymer  (polystyrene)  photonic  crystals.  Highly 
controlled  defect  structures  of  different  types  were  created  inside  three-dimensional 
polymeric  photonic  crystals  with  gaps  centered  in  the  visible  region.  Two  of  these 
structures  are  presented  in  Figure  .... 

The  main  purpose  of 
this  investigation 
was  to  determine  the 
optical  properties  of 
samples  containing 
intentional  defects. 


In  our 

experiments  we  used 
two  methods  to 
prepare  photonic 
crystals:  the  vertical 
deposition  method  and 
the  cell  method. 
Comparing  these  two 

types  of  colloidal  photonic  crystals  we  found  that  after  washing  away  the  unpolymerized  parts  of 
the  resin  subsequent  to  the  lithographic  process,  the  quality  and  the  stability  of  the  photonic 
crystal  made  using  the  cell  method  were  much  better.  Our  fabricated  structures  were  of  high 
enough  quality  to  guide  light  with  the  losses  during  guiding  of  the  fluorescence  emission  through 
the  defect  of  about  0.06dB/pm.  Better  matching  of  refractive  index  of  the  polymerized  resin  to 
that  of  the  polymer  spheres  forming  the  photonic  crystal  and  improvements  in  the  lithographic 
process  should  lead  to  lowering  the  losses. 


Figure  17.  A  illustrates  gratings  with  a  resolution  of  1pm, 
fabricated  using  two  photon  lithography,  10pm  below  the 
surface  of  the  crystal  made  using  the  vertical  deposition 
method.  Figure  2B  shows  a  beam  splitter  with  a  resolution 
of  1pm,  fabricated  3pm  below  the  surface  of  the  crystal. 
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IX.  Multibeam  Holographic  Lithography 

Collaboration:  P.  N.  Prasad,,  University  at  Buffalo; 

E.  Thomas,  MIT,  Boston 

Holographic  lithography  is  a  method  to  fabricate  periodic  structures  by  exposing  a 
photosensitive  precursor  to  a  defined  microstructure  obtained  by  interference  of  coherent 
light.  Utilizing  holographic  lithography  as  a  technique  for  the  fabrication  of  photonic 
crystals  has  some  advantages  over  the  other  proposed  methodologies.  For  instance, 
colloidal  crystals,  and  structures  made  by  using  colloidal  crystals  as  templates  are 
relatively  cheap  route  to  obtain  photonic  crystals.  However,  the  growth  time  is  relatively 
longer  and  the  structures  have  intrinsic  defects  such  as  stacking  faults,  vacancies  etc. 
More  importantly  they  are  only  suitable  for  the  fabrication  of  photonic  crystals  with  face- 
centered-cubic  (fee)  symmetry.  In  holographic  lithography,  however,  various  geometries 
with  sub-micron  resolution  can  be  defined  easily  and  written  on  a  photo-polymerizable 
resin,  which  makes  this  method  very  suitable  for  the  growth  of  photonic  band  gap 
materials. 

Our  team  (Prasad,  Thomas)  has  developed  a  procedure  for  fabricating  excellent  quality 
photonic  crystals,  in  high  refractive  index  photo-polmerizable  materials  by  using  four- 
beam  interference  lithography.  We  used  an  “umbrella”  configuration  with  two  different 
polarizations  of  the  recording  beams  to  fabricate  photonic  crystal  structures  in  a 
photosensitive  resins  In  the  umbrella  geometry,  three  side  beams  form  a  tetrahedron 
whose  base  is  parallel  to  xy  plane  in  local  xyz  Cartesian  coordinates.  The  fourth,  central 
beam,  is  the  tetrahedron’s  height  and  intersects  with  all  three  other  beams  1,  2  and  3  at 
the  tetrahedron’s  apex  point. 

Holographis  lithography  was  used  to  record  3D  photonic  crystal  structures  in  two 
different  resins:  ORMOCER  and  SU-8. 
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Figure  18  .Representative  SEM  images  of  photonic  crystal  structures 
produced  in  SU8  resin  using  the  umbrella  configuration. 


As  can  be  seen  from  figure  18  ,  crystalline  structures  with  a  very  long  range  order  can  be 
fabricated  by  holographic  lithography.  The  crystals  prepared  by  this  method  showed 
high  structural  qualities.  Combination  of  ability  to  produce  desirable  geometries  and  to 
fabricate  single  crystalline  structures  with  no  defects  makes  this  method  a  very  suitable 
route  for  the  fabrication  of  photonic  crystals. 


X.  Photonically  Aligned  Nanoarray 

Collaboration:  P.  N.  Prasad,  A.  N.  Cartwright,  University  at  Buffalo; 

R.  Vaia  and  T.  Running,  AFRL,  Dayton 

A  novel  promising  photonic  approach  for  nanofabrication  is  ordering  of 
nanoparticles  that  is  produced  by  light-driven  migration  of  nanoparticles.  This  method 
was  developed  at  AFRL,  Dayton  by  Bunning  and  Vaia.  Their  collaboration  with  us  has 
been  instrumental  in  utilizing  this  method  for  achieving  enhanced  performance  in  a 
variety  of  photonic  functions. 

This  method  utilizes  randomly  distributed  nanoparticles  in  a  matrix  which  can  be 
a  monomer  of  appropriate  structure  or  a  polymer.  In  a  holographic  approach,  when  two 
beams  are  crossed  to  generate  an  intensity  grating,  intensity  modulation  producing 
alternate  bright  and  dark  regions  occur.  A  photoprocess  induced  by  this  intensity 
modulation  pattern  could  be  started  with  a  monomer  which  photopolymerizes.  Where 
there  is  a  bright  light  it  forms  a  polymer  and  where  there  is  dark,  it  remains  as  a 
monomer.  One  can  cross-link  a  polymer  in  the  bright  light,  whereas  in  the  dark  region, 
no  cross-link  of  the  polymer  occurs.  If  the  nanoparticles  do  not  like  what  is  being 
formed,  they  will  be  squeezed  out  and  go  to  the  dark  region.  A  spatial  movement  of  the 
nanoparticles  will  occur  to  align  them  in  the  dark  regions,  producing  an  array  of 
nanoparticles. 

One  can  utilize  this  approach  to  align  an  array  of  Q-dots  and  metallic 
nanoparticles  or  to  produce  photonics  crystals  where  the  lined  nanodots  with  high 
refractive  index  produce  a  high  refractive  index  contrast.  Alternatively,  one  can  utilize 
the  polymer-dispersed  liquid  crystal  to  align  liquid  crystal  nanodroplets. 

In  collaboration  with  AFRL,  we  utilized  the  holographically  aligned  polymer 
dispersed  liquid-crystal  (H-PDLC)  grating  film  as  an  angle-dependent  and  narrow 
spectral-band  feedback  control  element  for  two-photon  pumped  lasing  in  a  dye  solution 
(APSS  in  DMSO).  The  grating  film  contained  an  ~1 5  pm  thick  polymer  film,  featuring  a 
maximum  reflectance  of  75%  at  561-nm  position  with  an  ~9-nm  spectral  bandwidth.  The 
output  lasing  wavelength  could  be  tuned  from  561.5  to  548.5  nm  and  the  lasing 
bandwidth  changed  from  5  to  3  nm  when  the  incidence  angle  on  the  grating  film  varied 
from  0°  to  22°.  The  overall  lasing  efficiency  was  measured  to  be  10%.  The  major 
advantages  of  using  a  H-PDLC  film  grating  are  (i)  the  possibility  of  integrating  the  lasing 
medium  with  the  cavity  feedback  and  spectral  tunable  element  together  to  form  a 
compact  waveguide  or  integrated  micro-lasing  device,  (ii)  the  lasing  tunability  by 
applying  a  external  electric  field. 

More  recently,  we  have  used  a  grating,  produced  holographically  by  this  method, 
as  a  one-dimensional  photonic  crystal  and  demonstrated  dramatic  enhancement  of  third- 
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harmonic  generation  in  cases  where  the  fundamental  wavelength  is  near  the  stop-gap  or 
when  the  third-harmonic  wavelength  is  near  the  stop-gap.  This  is  a  considerably 
significant  achievement  in  combining  both  these  features  to  achieve  practically  viable 
third-harmonic  generation. 


XI.  Nanostructured  Multicomponent  Erbium  Doped  Glass  for  Significantly 
Enhanced  Optical  Amplification 

Collaboration:  P.  N.  Prasad,  A.  N.  Cartwright  and  H.  Luo, 

University  at  Buffalo 

Non-DURINT  Team  Member:  Dr.  J.  O’Reilly,  University  at  Buffalo; 

T.  Smith,  Rochester  Institute  of  Technology 

Silica  glass  is  a  highly  attractive  matrix  for  optical  amplification  due  to  its 
favorable  thermal  and  mechanical  properties.  Rare  earth  doped  silica  glasses  are 
potential  materials  for  solid  state  fiber  lasers,  amplifiers,  and  upconversion  lasers.  At  the 
Institute  for  Lasers,  Photonics  and  Biophotonics,  a  process  has  been  developed  to  control 
and  modify  the  local  environment  of  the  doped  rare  earth  ions.  Nanostructure  control 
was  successfully  used  for  two  purposes:  1)  to  reduce  the  non-radiative  quenching  due  to 
hydroxyl  groups;  and  2)  to  allow  a  high  number  density  of  erbium  without  aggregation 
and  concentration  quenching.  These  advances  have  produced  an  enhanced  emission. 
Another  major  advancement  from  these  procedures  is  a  longer  lifetime  and  greatly 
enhanced  calculated  gain  at  1.5  pm.  A  provisional  patent  has  been  filed  for  this  process. 

More  recently,  we  have  concentrated  on  developing  Er  doped  nanostructured 
tellurite  glasses  which  have  high  refractive  index  and  lower  frequency  phonons  to  further 
enhance  amplification  properties  and,  at  the  same  time,  increase  band  width.  We  have 
used  a  novel  sol-gel  precursor  route  for  producing  tellurite  glasses  and  already  have  some 
initial  success  in  producing  films  of  tellurite  glasses. 


XII.  Optical  Microfabrication  of  Highly  Reflective  Volume  Bragg  Gratings  using 

Polymer  Nanocomposites 

A.N.  Cartwright,  P.N.  Prasad,  University  at  Buffalo 
T.  Burning,  R.  Vaia,  AFRL,  Dayton 

This  task  was  aimed  towards  developing  highly  efficient  and  simple  methods  of 
producing  reflection  grating  using  nanocomposites.  Recently,  in  collaboration  with  Tim 
Bunning  at  Wright  Patterson  Air  Force  Base,  we  demonstrated  that  organic  polymeric 
materials  can  yield  high  efficiency  reflection  gratings  through  a  fast  and  inexpensive 
holographic  photopolymerization  technique  Specifically,  multi-layer  grating  structures 
were  prepared  by  a  holographic  technique  using  photopolymer  syrups  (60wt  %) 
composed  of  monomer  (Dipentaerythritol  hydroxypentaacrylate),  photoinitiator  (Rose 
Bengal),  co-initiator  (N-phenylglycine),  cross-linking  monomer  (N-vinylpyrrolidinone), 
liquid  crystal  (TL213)  and  non-reactive  solvent  (toluene)  (20wt  %).  For  the  sample  with 
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silane  incorporated,  toluene  contains  1.5wt  %  3-aminopropyltriethoxysilane.  The  liquid 
crystal  is  commercially  available  from  Merck,  the  photoinitiator  was  obtained  from 
Spectra  Group  Limited,  and  all  other  chemicals  were  obtained  from  Aldrich.  The  pre¬ 
polymer  syrup  was  mixed  to  ensure  homogeneity  with  a  mixer  and  a  sonicator  and  then 
sandwiched  between  two  glass  slides.  Holography  was  performed  using  the  514  nm  line 
of  an  Ar  ion  laser  with  100  mW  of  power;  with  an  exposure  time  was  30  seconds.  After 
laser  exposure,  the  sandwiched  cell  was  post-cured  under  a  75  W  lamp  for  24  hrs.  The 
reflection  notch  could  be  observed  after  removing  the  cover  glass  slide  from  the  cell  and 
allowing  the  incorporated  solvent  to  evaporate.  Within  the  past  six  months,  we  have 
demonstrated  that  it  is  possible  to  create  porous  polymeric  gratings  without  using  liquid 
crystals.  Specifically,  we  have  demonstrated  that  a  mixture  of  monomer,  photoinitiator, 
co-initiator  and  formamide  provides  results  that  are  almost  identical  to  the  case  described 
in  detail  above  (formamide  replaces  the  liquid  crystal  and  acetone  combination).  This 
improved  technique  provides  a  much  simple  fabrication  method  for  producing  2D 
transmission  and  ID  reflection  gratings. 

In  addition,  we  have  demonstrated  that  chemical  vapors  can  easily  permeate  into 
the  voids  in  the  grating  film  and  consequently  change  the  average  RI  and  the 
corresponding  reflective  wavelength.  These  changes  in  optical  response  can  be  readily 
observed  due  to  the  large  index  difference  between  the  air  voids  (nair=T)  and  the  chemical 
vapor  (normally  >1.36).  Thus,  these  holographic  porous  polymer  gratings  can  act  as 
colorimetric  vapor  sensors.  The  wavelength  shifts  observed  in  these  structures  are  far  too 
large  to  be  explained  simply  by  the  displacement  of  air  in  the  pores  by  condensed 
acetone,  as  occurs  in  analogous  porous  silicon  vapor  sensors.  According  to  the  Bragg 
equation,  and  assuming  a  polymer  refractive  index  of  about  1 .5,  a  structure  with  200  nm 
grating  spacing  and  40%  porosity  would  exhibit  a  wavelength  shift  of  only  about  60  nm 
if  the  pores  were  completely  filled  with  liquid  acetone  (index  =  1.36).  This  suggests  that 
there  is  another  reversible  mechanism  occurring  that  affects  the  optical  thickness  nme A 

of  the  gratings.  One  possible  cause  for  this  could  be  absorption  of  acetone  solvent  by  the 
polymer,  which  would  cause  the  grating  to  swell,  increasing  the  physical  grating  spacing. 
Liquid  acetone  uptake  has  been  shown  to  substantially  swell  photochemically  crosslinked 
poly(n-butylacrylate),  increasing  the  diameter  of  the  polymer  network  by  about  80%. 

The  long-term  instability  of  the  porous  polymer  is  significantly  enhanced  by  adding  a 
small  amount  of  the  silane  compound,  3-aminopropyltriethoxysilane,  to  the  prepolymer 
syrup.  This  strongly  enhances  the  contact  between  the  polymer  and  the  glass,  while  not 
adversely  affecting  grating  formation.  In  contrast,  without  the  silane  compound  present, 
the  sample  would  easily  delaminate  from  the  glass  substrate,  when  immersed  in  the 
water.  It  should  be  noted  that  3-aminopropyltriethoxysilane  is  commonly  used  for 
preparing  amino  functional  surfaces  on  glass  and  silica.  This  suggests  that  the  use  of  this 
compound  could  aid  in  producing  porous  polymer  gratings  useful  for  biosensing 
applications,  both  by  rendering  such  devices  hydrophobic  and  by  providing  an  amino 
functional  group  within  the  porous  composite  which  could  potentially  be  used  to  target 
protein  analytes. 
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Optical  Micro  fabrication  of  Highly  Reflective  Volume  Bragg  Gratings  using 
Polymer  Nanocomposites 


In  collaboration  with  the  Air  Force  Research  Laboratory  at  Dayton,  a  novel 
approach  for  fabricating  wide  bandwidth  and  highly  reflective  Bragg  grating  structures 
(photonic  bandgap)  using  the  technique  of 
holographic  photopolymerization  of  liquid 
crystal  polymer  nanocomposites  was 
developed.  The  key  to  this  fabrication 
method  that  distinguishes  it  from  previous 
methods  is  the  use  of  a  non-reactive 
solvent,  acetone,  to  dissolve  the 
photoinitiator  and  coinitiator  in  an  acrylate 
monomer/liquid  crystal  (LC)  mixture.  The 
addition  of  acetone  results  in  the  creation 
of  controllable  periodic  voids  inside  the 
thin  film  after  the  acetone  evaporates.  The 
transmission  spectra  of  this  grating  together 
with  an  SEM  image  are  shown  in  Figure 
19.  Peak  reflectivity  as  high  as  80%  and  a 
broad  reflection  bandwidth  of  80  nm  were 
observed  in  the  reflection  gratings  formed 
with  acetone  present  in  the  starting 
mixture.  It  was  estimated  from  a  fit  to  the 
experimental  data  that  the  resulting  index 
mismatch  was  approximately  0.2  which  is  consistent  with  the  presence  of  air  voids.  This 
index  mismatch  is  very  high  when  compared  to  the  standard  photopolymer  system  for 
grating  formation.  By  backfilling  with  various  refractive  index  fluids,  tunable 
wavelength,  diffraction  efficiency,  and  bandwidth  reflection  notches  can  be  formed.  A 
provisional  patent  application  has  been  filed  for  this  technology. 
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Fig  19.  The  spectra  of  the  acetone 
sample  of  backfilling  with  different 
materials  (a)  air  voids  nair=l  (b) 
acetone,  na=1.36  (c)  chloroform 
nc=1.45  and  (d)  toluene  nt=1.50.  The 
inset  shows  SEM  of  the  grating 
structure  obtained  using  acetone  in  the 
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Prologue: 


Making  electro-optic  (EO)  modulators  and  other  EO  based  components  requires 
organizing  the  active  component  of  the  material  to  be  in  an  acentric  order.  The  individual 
components  should  have  as  large  an  EO  activity  and  be  as  well  ordered  as  possible  to 
make  the  best  device.  EO  activity  is  best  summarized  by  the  EO  coefficient,  r33 ,  which  is 


related  to  the  molecular  hyperpolarizability,  pzzz ,  the  number  density,  N ,  and  the  overall 


acentric  order  average  (or  parameter)  ^cos3  d'j  by  the  relation: 

r33  oc  Pz:zN  (cos3  6^j 


Therefore,  theoretical  understanding  of  the 
active  components  and  the  nature  of  the 
energetics  of  organization  of  those 
components  under  poling  fields  are  necessary 
to  design  devices  with  the  highest  EO  activity. 
Two  figures  of  merit  for  optical  activity  are 
plotted  as  a  function  of  time,  and  show  a 
Moore’s  Law  type  increase.  This  type  of 
progress  in  devices  has  been  accompanied  by 
theoretical  guidance  and  input  in  choosing 
optimal  components.  The  supporting 
theoretical  efforts,  which  have  aided  in 
designing  and  developing  novel  materials,  are 
detailed  below. 
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Goals  of  this  Durint: 


The  goal  of  this  work  has  been  to  fully  characterize  the  nature  of  organic-based  nonlinear 
optical  (ONLO)  materials.  This  consists  of  two  parts.  The  first  is  to  calculate  the 
hyperpolarizability  of  organic  chromophores,  which  are  the  individual  molecules  that  are 
the  active  constituents  of  ONLO  materials.  The  second  is  to  determine  the  net  order  of 
the  chromophores  when  incorporated  into  the  material  (as  a  condensed  matter  system).  In 
particular  the  computation  of  order  is  done  as  a  function  of  the  number  density  of  the 
chromophores.  The  work,  described  below,  has  succeeded  in  meeting  both  of  these 
goals. 

Supporting  Material  that  demonstrates  program  success: 

The  articles,  cited  below,  most  of  which  are  now  in  the  literature,  document  the  success 
of  the  effort  devoted  to  this  project.  One  very  important  result  has  been  to  compare 
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methods  of  computing  hyperpolarizability.  Many  literature  discussions  would  give  the 
impression  that  Density  Function  Theory  (DFT)  approaches  to  computation  do  not  give 
valid  results.  The  following  paper  shows  that  current  methods  using  either  DFT  or  FIF 
methods  will  give  accurate  estimates  of  molecular  hyperpolarizability: 

E.  R.  Davidson,  B.  E.  Eichinger,  B.  H.  Robinson,  “Hyperpolarizability:  Calibration  of 
Theoretical  Methods  for  Chloroform,  Water,  Acetonitrile,  and  p-Nitroaniline”;  Optical 
Materials,  In  Press,  Corrected  Proof,  Available  online  8  May  2006 

Fundamentally,  these  calculations  demonstrate  that  the  current  computational  methods 
are  indeed  able  to  quantitatively  and  accurately  (to  within  experimental  error)  predict  the 
hyperpolarizability  of  ONLO  chromophores,  when  in  the  gas  phase. 

The  following  article  makes  a  comparison  among  five  different  established  quantum 
mechanical  codes  to  determine  the  relative  accuracy  of  the  methods.  A  test  set  of  16 
ONLO  chromophores  were  studied.  The  relative  rank  order  of  the  results  was 
independent  of  method.  Some  disagreements  in  hyperpolarizability  were  found  among 
the  methods  but  they  were  considered  to  be  rather  modest.  The  conclusion  of  the  study  is 
that  one  may  choose  any  standard  method  to  compute  the  complete  molecular 
hyperpolarizability  tensor  of  individual  molecules  and  the  relative  results  will  be 
consistent  with  other  methods. 

C.  Isborn,  A.  LeClercq,  F.  D.  Vila,  J.  L.  Bredas,  B.  E.  Eichinger  B.  H.  Robinson, 
“Comparison  of  Hyperpolarizabilities  Calculated  with  Various  Quantum  Mechanical 
Methods”.  Under  Review  J.  Phys.  Chem.  A. 

Because  this  article  is  still  under  review, 
one  of  the  most  important  figures  in  it  is 
reproduced  here:  The 
hyperpolarizabilities  are  computed  for 
16  different  chromophores  of  interest  to 
experimentalists.  The  plot  is  a 
correlation  plot  for  the  different 
molecules.  Five  different  techniques 
were  used,  including  one  Hartree-Fock 
(HF)  method  and  one  Density  Functional 
Theory  (DFT)  method  (from  the 
Gaussian  package),  one  method  from 
the  Goddard  group,  one  DFT  method 
from  Accelrys,  and  one  semi-empirical  Intermediate  Neglect  of  Differential  Overlap 
(INDO)  method  from  the  laboratory  of  J.L.Bredas.  As  a  correlation  plot,  the  HF  method 
is  taken  as  the  standard  and  plotted  along  the  diagonal.  The  results  from  the  different 
methods  are  plotted  on  the  ordinate  for  the  corresponding  value  determined  by  HF 
method  along  the  abscissa.  .  The  hyperpolarizability  spans  three  decades  in  values.  The 
general  conclusion  is  that  all  methods  follow  the  diagonal  quite  well;  any  method,  if  used 
consistently  will  give  relative  results  that  can  be  used  to  predict  better  performance. 
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A  collaboration  between  the  Robinson  and  Rehr  groups  at  the  University  of  Washington 
has  resulted  in  a  paper  that  documents  a  new  method  for  computing  optical  properties  of 
ONLO  chromophores  with  improved 
accuracy: 

“Real  time  approach  for  non-linear 
optical  response  of  photonic 
molecules”,  Y.  Takimoto,  J.  J.  Rehr,  1 
F.  D.  Vila,  and  B.  E.  Eichinger,  1 
Phys.  Rev.  B  2006,  in  preparation. 

x 

The  method  uses  Time  Domain 
DFT.  This  results  in  more  accurate 
estimators  for  optical  responses.  The 
absorption  spectra  shown  here  for  a 
typical  chromophore  (FTC), 

compared  with  the  experimental  spectrum,  are  computed  as  Fourier  Transforms  of  the 
appropriate  time  domain  optical  function.  The  work  was  done  for  several  rotameric 
forms  of  FTC.  Quantum  mechanical  calculations  of  the  effect  of  solvents  and  varied 
rotameric  states  of  the  ONLO  chromophores  have  been  carried  out  and  are  documented 
in  the  following  articles: 

Kinnibrugh,  T.;  Bhattacharjee,  S.;  Sullivan,  P.;  Isborn,  C.;  Robinson,  B.  H.;  Eichinger,  B. 
E.,  “Influence  of  Isomerization  on  Nonlinear  Optical  Properties  of  Molecules”, 

J.  Phys.  Chem.  B.;  (Article);  2006;  110(27);  13512-13522. 

Liao,  Y.;  Eichinger,  B.  E.;  Firestone,  K.  A.;  Haller,  M.;  Luo,  J.;  Kaminsky,  W.;  Benedict, 
J.  B.;  Reid,  P.  J.;  Jen,  A.  K.-Y.;  Dalton,  L.  R.;  Robinson,  B.  H.  “Systematic  Study  of  the 
Structure-Property  Relationship  of  a  Series  of  Ferrocenyl  Nonlinear  Optical 
Chromophores”.  J.  Am.  Chem.  Soc.;  (Article);  2005;  127(8);  2758-2766. 

Liao,  Y.;  Bhattacharjee,  S.;  Firestone,  K.  A.;  Eichinger,  B.  E.;  Paranji,  R.;  Anderson,  C. 
A.;  Robinson,  B.  H.;  Reid,  P.  J.;  Dalton,  L.  R.,  “Antiparallel-Aligned  Neutral-Ground- 
State  and  Zwitterionic  Chromophores  as  a  Nonlinear  Optical  Material”,  J.  Am.  Chem. 
Soc.;  (Article);  2006;  128(21);  6847-6853. 

In  addition  to  the  optical  spectrum  of  FTC,  other 
optical  properties  can  also  be  computed  as  a 
function  of  the  wavelength  of  the  light  field. 

Here  is  the  hyperpolarizability  (along  the 
molecular  z  axis)  as  a  function  of  the 
wavelength  of  light.  The  vertical  lines  represent 
two  important  telecom  wavelengths.  Notice  that 
the  hyperpolarizability  nearly  vanishes  at  one  of 
the  telecom  wavelengths  and  it  close  to  the  long 
wavelength  limit  at  the  other  telecom 
wavelength.  This  underscores  the  need  to 
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understand  the  wavelength  dependence  of  the  full  hyperpolarizability  tensor. 

Marks  and  Ratner  have  suggested  that  a  point  of  instability  in  a  chromophore  can  lead  to 
exceedingly  large  molecular  hyperpolarizability.  They  used  a  Multi-Reference  (MR) 
INDO  type  method  to  show  that  a  twisting  about  a  double  bond  at  nearly  90°  may  lead  to 
considerable  increase  in  hyperpolarizability.  A  concept  of  this  potential  import  deserves 
to  be  treated  theoretically  with  the  best  possible  methods.  Therefore,  we  undertook  to 
determine  whether  more  exacting  quantum  mechanical  calculations  would  support  this 
conclusion.  We  used  HF  (CSS,  MP2  and  MP4)  and  DFT  (several  functionals  and 
broken  spin)  methods  and  found  almost  no  evidence  to  support  the  original  claim  of 
Marks  and  Ratner.  We  hope  to  be  proven  incorrect  by  experiment. 

Christine  M.  Isborn,  Ernest  R.  Davidson,  Bruce  H.  Robinson,  “Ab  Initio 
Dirad ical/Zwitterionic  Polarizabilities  and  Hyperpolarizabilities  in  Twisted  Double 
Bonds”.  J.  Phys.  Chem.  A.;  (Article);  2006;  1 10(22);  7189-7196 

More  recently  we  have  examined  whether  or  not  different  donor  and  acceptor  groups 
would  improve  hyperpolarizability.  Below  we  illustrate  that  a  potential  40  fold 
improvement  over  the  currently  used  chromophore,  FTC,  in  hyperpolarizability  is 
possible  theoretically,  with  some  red-shift  in  the  absorbance  band.  Working  with 
experimentalists,  we  are  now  synthesizing  several  of  these  compounds.  This  work  has 
not  been  published  as  it  is  in  the  preliminary  stages  and  remains  to  be  seen  whether  DFT 
has  any  predictive  power  here. 
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The  second  major  thrust  has  been  to  understand  the  energetics  of  the  organization  of  the 
chromophores  in  condensed  phases.  The  need  for  this  is  illustrated  and  underscored  by 
the  now  well  known  observation  that  the  EO  response  (or  activity)  depends  very  strongly 
on  the  number  density  of  chromophores,  when  the  chromophores  are  contained  in  a  host 
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material.  The  figure  (to  the  right)  shows 
that  the  EO  activity  increases  with 
chromophore  number  density  but  then 
drops  off.  We  attribute  the  fall-off  to 
intermolecular  electrostatics.  The  large 
dipolar  interactions  between 
chromophores  tend  to  generate  an 
interaction  that  mitigates  against 
alignment  with  the  poling  field.  .  The  fit 
to  the  data  in  this  figure  is  done  with  a 
model  that  takes  into  account  the 
intermolecular  electrostatic  interactions  as 
well  as  the  interactions  of  the  individual 
chromophores  with  the  poling  field,  as  well  as  the  prolate  shape  of  the  chromophores. 
The  model  considers  only  the  simplest  phenomenological  theory.  However,  many  of  the 
features  of  the  model  have  been  verified  and  amplified  on  by  detailed  Monte  Carlo 
simulations  of  systems  containing  many  ellipsoids.  The  following  paper  documents  this 
effect  for  chromophores  represented  by  simple  ellipsoids  with  different  shapes. 

Nielsen,  R.  D.,  Rommel,  H.  L.  and  Robinson,  B.  H.  "Simulation  of  the  Loading 
Parameter  in  Organic  Nonlinear  Optical  Materials".  Journal  of  Physical  Chemistry  B 
2004;  108:  8659-8667. 
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The  following  discussion  about  the  order  in  condensed  matter  systems  simulated  by 
Monte  Carlo  methods  has  not  been  published.  Therefore  a  rather  lengthy  discussion  of 
the  newest  results  will  be  given.  We  have  continued  to  examine  the  effects  of  the 
intermolecular  potential  on  the  total  acentric 

ordering  that  is  generated  by  the  presence  of  . . . . . ' . ^ 

a  poling  field.  In  particular  we  now  carry  J  * '  0',“" 

out  off-lattice  calculations  (which  means  the 
ellipsoids  are  allowed  to  both  translate  and 
rotate)  and  NPT  as  well  as  NVT  calculations 
that  determine  the  natural  density  of  the 
chromophore  system.  We  have  found  a 
surprising  result:  When  the  ellipsoids  are 
allowed  to  move  off  lattice  the  effect  of  the 
shape  of  the  ellipsoids  on  the  order  is  nearly 
nonexistent.  The  figure  to  the  right 
illustrates  that  for  spherical,  oblate  and 
prolate  ellipsoids  with  the  same  volume,  the 
loading  parameter  (which  is  the  product  of  the 
order  parameter  and  the  number  density)  are  all 
very  similar. 
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advances  and  improvements  on  the  theory.  In  particular  we  have  modeled  the  detailed 
molecular  structure  of  the  molecules  by  a  set  of  ellipsoids  that  map  rotations  of  ellipsoids 
directly  to  bonds.  Therefore,  after  moving  the  ellipsoids  about  the  connections  between 
them  (where  the  small  green  circles  are  on  the  molecules) ,  it  is  possible  to  reconstruct 
the  molecule.  Moreover,  the  extrapolation  of  molecules  to  ellipsoids  does  not  allow 
molecularly  rigid  entities  to  rotate.  To  illustrate  how  the  ellipsoids  are  constructed  we 
show  a  simple  chromophore-containing  molecule  with  the  ellipsoids  overlaid  on  it.  The 
chromophore  part  (that  is  molecularly  rigid)  is  represented  by  two  ellipsoids  that  cannot 
move  with  respect  to  each  other.  However,  CH2  moieties  are  able  to  move  with  respect 
to  their  neighbors  but  only  about  the  bonds  connecting  them.  In  this  way  the  complex 
nature  of  the  molecule  is  reduced  to  a  simplified  united-atom  or  coarse  grained  approach. 


The  Dalton  group  has  made  several  tri-chromophore  containing  molecules  (synthesized 
using  the  principles  of  dendritic  synthesis).  We  have  recently  modeled  the  EO  activity 
and  now  compare  the  theoretical  results  with  those  of  experiment.  The  figure  shows  the 
structure  of  the  system  containing  over  thirty  tri-chromophric  molecules  at  the  density 
expected  for  this  system  (which  is  around  1  gram/cc).  There  is  order  in  this  system  due 
to  the  presence  of  the  poling  field.  In  the  absence  of  a  poling  field  order  is  never 


generated. 

The  rainbow  colored  parts  of  the  molecules  are  the  chromophores.  The  blue  end 
represents  the  donor  end  of  the  chromophore.  There  is  net  order  (around  0.25)  due  to  the 
poling  field.  At  this  density  the  chromophore  interactions  (through  the  inter- 
chromophore  dipole-dipole  interactions) 

reduce  the  overall  order.  1  . . .  "  1 


The  order  of  these  tri-chromophore 
containing  systems  is  compared  to  that  of 
simple  (off  lattice)  spheres  at  comparable 
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Chromophore  Number  Density  (N/cc) 


chromophore  density.  One  can  see  that  even  at  low  density  the  tri-chromophore 
molecules  have  low  order.  That  order  does  not  change  (substantially)  with  density.  This 
happens  because  three  chromophores  are  always  in  close  proximity  even  at  low  density. 
Both  molecules  (PAS41  and  PAS33),  using  the  same  root  chromophore,  have  very 
similar  order.  PAS33  has  a  slightly  higher  chromophore  density  (at  1  gram/cc)  because 
the  molecule  contains  less  additional  non-chromophore  material.  The  theory  predicts  that 
the  more  densely  packed  molecule  will  have  about  twice  the  EO  activity.  The 
experimental  results,  shown  below,  are  that  the  more  densely  packed  tri-chromophore 
system  is  about  50%  more  active.  This  is  surprising  because  an  increase  in  molecular 
density  usually  causes  an  increase  in  intermolecular  interactions  which  reduces  the 
loading.  Here  however  the  loading  is  increased  by  about  80%  theoretically  and  50% 
experimentally.  The  details  of  the  molecular  architecture  appear  to  be  important  causing 
this  reversal  in  the  general  trend.  The  error  bars  shown  in  the  figure  are  not  experimental 
error  as  these  are  theoretical  results.  They  represent  the  range  of  values  of  the  order 
parameters  encountered  amongst  all  of  the  replicate  runs  (about  50,000  such  runs  at  each 
density).  The  individual  Monte  Carlo  runs  yield  individual  results  because  the 
organization  is  determined  by  different  random  (technically,  pseudo-random)  numbers. 

To  the  right,  the  EO  activity  of  the  fundamental  chromophore  in  APC  (shown  as  Xs)  is 
plotted  as  a  function  of  the  density  of 
chromophores.  The  fit  of  the  model  (solid 
squares)  to  the  data  has  no  adjustable 
parameters.  The  hyperpolarizability  has 
been  measured  by  HRS.  The  agreement 
with  the  system  using  the  polymer,  APC,  as 
a  host  is  excellent.  The  agreement  with  the 
tri-chromophore  systems  (Os)  is  not  as 
good.  As  shown,  the  general  trend  of  the 
two  tri-chromophore  molecules  is  about 
50%  larger  than  that  predicted  by  the 
spherical  ellipsoids  and  about  twice  that 
predicted  using  the  tri-chromophore 
molecules  (see  the  previous  figure).  The 
fact  that  our  results  are  within  a  factor  of  two  of  the  experimental  system  is  really  quite 
remarkable.  The  details  of  the  intermolecular  interactions  have  been  replaced  by  very 
simple  rules  about  the  van  der  Waals  interactions  and  simple  dipoles  at  the  center  of  the 
chromophores.  It  is  possible  that  the  actual  dipole  moments  and  hyperpolarizability  are 
somewhat  different  in  the  systems  that  are  purely  tri-chromophore  molecules  as  opposed 
to  the  case  where  the  individual  chromophores  are  doped  into  APC.  Overall,  the  ability 
to  simulate  large  complicated  systems  is  really  very  gratifying.  By  contrast,  one  cannot 
obtain  similar  results  using  atomistic  force  fields  and  commercially  available  codes:  the 
forces  are  so  large  in  those  systems  and  codes  that  the  chromophores  do  not  align  or 
move  under  experimental  conditions. 

The  codes  to  do  the  Monte  Carlo  simulations  are  written  in  Matlab  and  are  available  upon 
request. 
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Papers  that  cite  the  support  of  AFQSR  (some  are  included  in  the  above  discussion): 


Liao,  Y.;  Firestone,  K.  A.;  Bhattacharjee,  S.;  Luo,  J.;  Haller,  M.;  Hau,  S.;  Anderson,  C. 
A.;  Lao,  D.;  Eichinger,  B.  E.;  Robinson,  B.  H.;  Reid,  P.  J.;  Jen,  A.  K.-Y.;  Dalton,  L.  R. 
“Linear  and  Nonlinear  Optical  Properties  of  a  Macrocyclic  Trichromophore  Bundle  with 
Parallel-Aligned  Dipole  Moments”.  J.  Phys.  Chem.  B.;  (Article);  2006;  1 10(1 1);  5434- 
5438. 

Nielsen,  R.  D.,  Rommel,  H.  L.  and  Robinson,  B.  H.  "Simulation  of  the  Loading 
Parameter  in  Organic  Nonlinear  Optical  Materials".  Journal  of  Physical  Chemistry  B 
2004;  108: 8659-8667. 

Christine  M.  Isborn,  Ernest  R.  Davidson,  Bruce  H.  Robinson,  “Ab  Initio 
Diradical/Zwitterionic  Polarizabilities  and  Hyperpolarizabilities  in  Twisted  Double 
Bonds”.  J.  Phys.  Chem.  A.;  (Article);  2006;  110(22);  7189-7196 

Kinnibrugh,  T.;  Bhattacharjee,  S.;  Sullivan,  P.;  Isborn,  C.;  Robinson,  B.  H.;  Eichinger,  B. 
E.,  “Influence  of  Isomerization  on  Nonlinear  Optical  Properties  of  Molecules”, 

J.  Phys.  Chem.  B.;  (Article);  2006;  110(27);  13512-13522. 

Liao,  Y.;  Eichinger,  B.  E.;  Firestone,  K.  A.;  Haller,  M.;  Luo,  J.;  Kaminsky,  W.;  Benedict, 
J.  B.;  Reid,  P.  J.;  Jen,  A.  K.-Y.;  Dalton,  L.  R.;  Robinson,  B.  H.  “Systematic  Study  of  the 
Structure-Property  Relationship  of  a  Series  of  Ferrocenyl  Nonlinear  Optical 
Chromophores”.  J.  Am.  Chem.  Soc.;  (Article);  2005;  127(8);  2758-2766. 

Liao,  Y.;  Bhattacharjee,  S.;  Firestone,  K.  A.;  Eichinger,  B.  E.;  Paranji,  R.;  Anderson,  C. 
A.;  Robinson,  B.  H.;  Reid,  P.  J.;  Dalton,  L.  R.,  “Antiparallel-Aligned  Neutral-Ground- 
State  and  Zwitterionic  Chromophores  as  a  Nonlinear  Optical  Material”,  J.  Am.  Chem. 
Soc.;  (Article);  2006;  128(21);  6847-6853. 

“Real  time  approach  for  non-linear  optical  response  of  photonic  molecules”,  Y. 

Takimoto,  J.  J.  Rehr,  F.  D.  Vila,  and  B.  E.  Eichinger,  Phys.  Rev.  B  2006,  in  preparation. 

C.  Isborn,  A.  LeClercq,  F.  D.  Vila,  J.  L.  Bredas,  B.  E.  Eichinger  B.  H.  Robinson, 
“Comparison  of  Hyperpolarizabilities  Calculated  with  Various  Quantum  Mechanical 
Methods”.  Under  Review  J.  Phys.  Chem.  A. 

E.  R.  Davidson,  B.  E.  Eichinger,  B.  H.  Robinson,  “Hyperpolarizability:  Calibration  of 
Theoretical  Methods  for  Chloroform,  Water,  Acetonitrile,  and  p-Nitroaniline”;  Optical 
Materials,  In  Press,  Corrected  Proof,  Available  online  8  May  2006 
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To:  Prof.  Paras  PRASAD,  University  at  Buffalo 
Subject:  DURINT  Final  Report 

From:  Prof  Jean  Frechet  University  of  California,  Berkeley 
A.  Technical  Report  Information 

Working  in  close  collaboration  with  the  Wright  Patterson  team  led  by  Dr.  Loon-Seng 
Tan,  of  the  Polymer  Branch,  AFRL/MLBP,  Materials  &  Manufacturing  Directorate,  Air 
Force  Research  Laboratory  and  the  Buffalo  group  led  by  Professor  Paras  Prasad,  we  have 
designed  and  prepared  a  variety  of  photoactive  compounds  including  novel  multiphoton 
absorbing  chromophores  systems  based  on  a  dendritic  framework  as  well  as  micellar 
laser  responsive  systems.  The  multi-photon  light  harvesting  systems  can  harvest  and 
reprocess  the  light  energy  received  and  transform  it  in  another  form  of  energy  making 
them  useful  in  a  variety  of  technological  applications  and  medical  applications. 

Therefore,  the  energy  absorbed  by  the  multiple  2-photon  chromophores  present  in  the 
system  can  be  focused  onto  a  single  “receiver”  core  where  it  is  reprocessed.  We  have  also 
studied  the  fundamental  photophysical  behavior  of  these  systems  as  well  as  their  uses  in 
applications  such  as  up-conversion  of  energy,  two-photon  pumped  lasing,  optical  power 
limiting,  photocatalytic  deactivation  of  chemical  agents,  targeted  photodynamic  therapy, 
and  laser  induced  release  of  active  compounds.  Other  potential  applications  not  studied 
under  this  program  due  to  time  and  funding  limitations  involve  sensing  and  field-use 
power  sources.  In  order  to  make  the  systems  we  developed  practical  for  use  with 
inexpensive  methods  of  application  we  have  made  our  multiphoton  antennae  and  related 
functional  and  light  triggered  systems  soluble  in  common  organic  solvents,  thus  enabling 
their  use  with  standard  spin-coating  equipment. 

For  example  we  have  prepared  a  novel  fully  modular  system  consisting  of  a 
energy-harvesting  system  containing  a  porphyrin  (or  phthalocyanin)  moiety,  2-photon 
absorbers  and  solubilizing  moieties  for  use  in  optical  power  limiting  or  in  the  generation 
of  highly  reactive  singlet  oxygen.  This  modular  system  is  designed  to  be  powered  by  IR 
laser  excitation  of  peripheral  two-photon  chromophores;  it  also  includes  solubilizing 
groups  since  both  the  core  porphyrin  and  the  two-photon  chromophores  tend  to  exhibit 
very  low  solubilities.  We  have  explored  metallation  of  the  central  porphyrin  to  tune  the 
properties  of  this  system  and  optimize  its  photophysical  properties  depending  on  the 
metal.  For  example,  the  porphyrin-silver  complex  performs  well  in  optical  power 
limiting  applications,  as  it  is  capable  of  converting  high  intensity  red  laser  light  to  heat 
via  non-radiative  decay.  A  schematic  representation  of  a  porphyrin  centered  optical 
power  limiting  dendrimer  is  shown  in  Figure  1.  In  contrast  if  palladium  is  used  in  the 
porphyrin  core,  the  resulting  light  harvesting  complex  becomes  a  remarkable 
photocatalyst  capable  of  near  100%  efficient  production  of  singlet  oxygen  upon  two- 
photon  irradiation.  A  water-soluble  version  of  this  system  with  AF608-341  two-photon 
chromophores  was  designed  and  synthesized  in  order  to  also  demonstrate  two-photon 
excited  FRET  in  aqueous  medium;  such  a  system  is  shown  in  Figure  2.  Extending 
nonlinear  optical  phenomena  to  aqueous  media  can  broaden  the  scope  of  our  modular 
energy  harvesting  antennae  to  biologically  or  field-relevant  applications  and  other 
aqueous-based  applications. 

We  have  also  designed  novel  materials  for  two-photon  pump  lasing,  but  much  of 
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this  work  had  to  be  discontinued  to  the  limited  manpower  available  and  a  focus  on  newly 
designed  optical  power  limiting  materials. 

Laser  technology  has  resulted  in  a  myriad  of  applications  both  in  the  civilian  and 
military  markets.  The  development  of  a  portable  laser  system  capable  of  inflicting  severe 
retinal  and  electronic  sensor  damage  for  use  in  military  applications  has  been  reported. 
Therefore  there  is  an  increased  need  for  protection  from  such  high  energy,  frequency 
agile  lasers  to  even  simple  laser  systems  such  as  hand-held  laser  pointers.  Optical  power 
limiting  materials  (OPL)  offer  the  possibility  of  maintaining  high  transparency  at  ambient 
light  intensities  while  maintaining  efficient  and  instantaneous  protection  against  the 
intense  energy  delivered  by  high  power  lasers. 

Exploring  a  new  design  in  a  proof-of-concept  study,  we  have  synthesized  a  novel 
OPL  organic  material  1  (Figure  3),  which  contains  a  two-photon  absorbing  chromophore 
(TPAC)  and  a  charge  separation  (CS)  component  comprised  of  perylene  tethered  to 
poly(3-hexyl  thiophene)  (P3HT).  The  design  of  this  system  allows  for  fluorescence 
resonance  energy  transfer  (FRET)  from  the  TPAC  to  the  CS  component  resulting  in  a 
charge-separated  state  localized  on  perylene  and  P3HT.  It  is  known  that  when  perylene 
and  P3HT  derivatives  are  covalently  bound  they  undergo  rapid  (~2  ps)  photoinduced 
electron  transfer  resulting  in  a  charge  separated  species  when  either  chromophore  is 
excited,  followed  by  a  charge  recombination  that  is  more  than  an  order  of  magnitude 
slower  (~50  ps).  Evidence  for  this  process  can  be  observed  in  the  quenched  fluorescence 
emission  of  both  chromophores.  Transient  absorption  of  this  charge-separated  species 
occurs  within  the  NIR  region.  In  order  to  enhance  the  OPL  ability  of  this  charge 
separation  dimer  for  two-photon  excitation  within  the  NIR  spectral  region,  we  utilized  the 
concept  of  two-photon  excited  FRET  from  a  known  TPAC  to  access  this  transient  charge 
separated  species.  In  this  way  we  can  harness  the  two-photon  absorption  efficiency  of 
the  AF-341  derivative  and  funnel  its  excited  state  energy  into  producing  the  transient 
charge  separated  species  required  for  reverse  saturable  absorption  (RSA)  and  ultimately 
OPL  of  the  input  laser  pulse.  Both  charge  separation  components  and  two-photon 
absorbing  chromophores  have  been  utilized  independently  for  optical  power  limiting 
(OPL),  but  we  have  found  that  combining  these  in  a  single  macromolecule  provides 
enhanced  OPL  behavior  relative  to  the  individual  components.  This  work  has 
demonstrated  the  synthesis  of  a  novel  optical  power  limiting  material  containing  O- 
alkylated  AF-341  as  the  two-photon  absorbing  chromophore  covalently  appended  to  a 
perylene-P3HT  charge-separation  dimer.  This  system  is  transparent  within  the  NIR 
region  under  low  power  intensities  and  then  becomes  strongly  absorbing  at  high  laser 
intensities.  Such  optical  properties  may  offer  protection  of  NIR  sensors/heat  sensors 
against  damage  from  high  intensity  light.  One-  and  two-photon  excitation  studies  indicate 
that  FRET  occurs  from  the  TPAC  to  the  perylene-P3HT  dimer  resulting  in  charge 
separation.  To  our  knowledge,  this  is  the  first  reported  organic  material  capable  of  TP 
excited  FRET  to  a  CS  dimer.  It  is  believed  that  a  mechanism  of  absorption  derived  from  a 
rapidly  photogenerated  charge  separated  state  is  what  gives  rise  to  the  excellent  observed 
optical  power  limiting  behavior. 

Similar  approach  could  be  used  to  design  novel  multichromophoric  light 
harvesting  molecules  that  absorb  and  convert  UV  and  visible  radiation. 

During  the  course  of  this  study,  numerous  samples  of  active  materials  have  been 
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exchanged  with  Dr.  Loon  Seng  Tan  of  the  Air  Force  Research  Laboratories  at  Wright 
Patterson  AFB  and  numerous  measurements  and  model  applications  have  been  carried 
out  in  his  laboratory. 

Finally,  we  have  carried  out  studies  that  have  unambiguously  demonstrated  the 
feasibility  of  stimuli-responsive  micellar  systems  capable  of  releasing  an  active  payload 
upon  IR  laser  excitation.  This  is  normally  very  difficult  to  achieve  since  the  energy  of  IR 
photons  is  generally  much  too  low  to  enable  chemical  transformations.  However,  the 
combination  of  multiphoton  chromophores  than  can  “pool”  energy  to  achieve  sufficient 
energy  levels  (e.g.  pool  the  energy  of  2  or  more  IR  photons  to  achieve  an  energy  level 
equivalent  to  that  of  a  UV  photon)  with  a  moiety  capable  of  undergoing  a  drastic 
structural  change  upon  irradiation  has  allowed  us  to  demonstrate  the  feasibility  of  this 
approach  with  a  family  of  diazo-l,2-napthoquinones.  Micelles  produced  from 
appropriately  designed  polymeric  versions  of  this  2-photon  absorbing  chromophore  have 
been  shown  to  be  capable  of  a  photoinduced  molecular  rearrangement,  which  drastically 
alter  their  physical  properties,  while  enabling  the  IR  triggered  release  of  a  payload 
(Figure  4).  This  concept  of  multiphoton  induced  rearrangement  and  payload  release 
shows  promise  for  light  triggered  stimuli-responsive  systems  for  field  use. 

B.  Provide  program  statistics  as  outlined  below: 


(1)  Number  of  PI  and  Co-PI  involved  in  the  research  project : _ 1 

(list  on  an  attached  page) 

(2)  Number  of  Post  Doc  Supported  under  AFOSR: 

_ 3_ 

(list  on  an  attached  page) 

(3)  Number  of  graduate  students  supported  by  AFOSR:  _ 6 

(list  on  an  attached  page) 

(4)  Other  researchers  supported  by  AFOSR: 


(list  on  an  attached  page) 

(5)  Number  of  publications  by  Pi's  in  the  last  3  year  period  (2003-06)  in 

refereed  journals  _ 110 _ 
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defined  as  some  research  results  that  were  used  by  technology  developers 
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Figure  1.  AF608-341  and  metallated  derivatives  that  demonstrate 
efficient  energy  transfer  and  enhanced  photophysical  properties  for 

Optical  Power  Limiting. 
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Figure  2.  Water-Soluble  system  capable  of  multi-photon  excited  FRET 
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Figure  3.  Optical  power  limiting  via  1  or  2-photon  process. 
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Figure  4.  Laser  multiphoton  induced  degradation  of  a  robust  polymeric 
micelle  with  release  of  a  active  compound  (shown  as  green  dots). 
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Objectives 

The  objective  of  this  effort  was  to  explore  electronic  transport  in  molecular  systems. 

Status 


We  elucidated  the  major  transport  mechanisms  in  these  structures,  and  developed  a  new 
spectroscopic  technique  for  analyzing  these  devices.  We  also  discovered  a  new  sub¬ 
wavelength  optical  effect  in  metallic  nanowires  as  a  result  of  the  collaboration  through 
this  grant.  The  progress  is  detailed  in  the  next  section 

Accomplishments/New  Findings 

Background 

In  recent  years,  investigators  have  made  tremendous  strides  in  demonstrating  the 
basic  capabilities  of  molecular-based  electronics.  However,  this  work  consists  largely  of 
simple  circuits,  fabricated  under  carefully  controlled  laboratory  conditions,  by  extremely 
well  trained  investigators.  The  goal  of  this  program  is  to  develop  a  hybrid 
molecular/solid-state  circuit  configuration  that  will  be  fabricated  by  a  more  scaleable  and 
production-worthy  vapor  phase  fabrication  process.  The  hybrid  approach  builds  on  and 
combines  existing  capabilities  of  molecular  and  solid-state  circuits:  The  functional  unit  is 
a  bistable  molecular  gate  consisting  of  either  a  stacked  negative  differential  resistive 
(NDR)  or  memory  molecular  device.  This  device  has  already  been  successfully 
demonstrated.  These  active  molecular  gates  will  then  be  mated  with  a  MOS-based  post¬ 
processing  circuit  to  make  dense,  high-performance  cells.  This  union  will  be  facilitated 
by  new  vapor  phase  deposition  process  will  allow  vastly  more  controlled  fabrication  of 
the  fully  stacked  molecular  gates  upon  a  MOS-based  substrate  template. 

Current  work  on  molecular  electronics  relies,  almost  entirely,  on  liquid  phase 
processing.  Liquid  phase  processing  is  simple,  it  is  inexpensive,  and  it  builds  upon  a 
large  body  of  experience  with  self-assembled  molecular  monolayers  (SAMs).  We  now 
want  to  drive  molecular  circuits  to  a  complexity  that  will  match  (and  eventually  exceed) 
that  of  solid-state  microelectronics.  Early  microelectronic  processing  was  also  based  on 
the  use  of  liquids.  However,  as  circuit  complexity  and  area  increased,  liquid-based 
processing  could  not  deliver  the  uniformity  and  reproducibility  required  for  successful 
fabrication.  Thus,  virtually  all  modern  semiconductor  tools  are  based  on  the  use  of  gases. 
Gas  phase  processing  offers  far  superior  uniformity  and  reproducibility.  Indeed,  when 
surface  reactions  dominate,  it  is  possible  to  simultaneously  deposit  or  etch  the  surfaces  of 
dozens  of  stacked  wafers.  More  subtly,  gases  (and  ions)  allow  for  very  precise  and 
shallow  etching.  This  is  essential  in  the  fabrication  of  a  circuit  of  microprocessor-level 
complexity:  microprocessors  cannot  be  fabricated  in  a  single  deposition  step.  The 
operation  of  the  device  cell,  and  interconnection  of  cells,  inevitably  requires  a  complex 
multi-level  structure.  It  is  thus  essential  that  the  technology  allow  for  repeated  reinsertion 
into  deposition  systems.  For  such  a  sequence  of  depositions  to  work,  one  must  be  able  to 
clean  and  prepare  underlying  structures,  in-situ,  for  the  deposition  of  new  materials. 
Controlled  deposition  and  etching  are  thus  essential. 

This  project  is  a  collaboration  between  University  of  Virginia  and  Yale  University. 
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Yale’s  responsibility  was  the  fabrication  and  characterization  of  molecular  devices;  UVA 
had  the  responsibility  of  vapor  phase  molecular  deposition.  This  report  focuses  on  the 
fabrication  and  transport  data. 


Progress  and  Accomplishments 

The  suggestion1  and  demonstration2  of  utilizing  molecules  as  the  active  region  of 
electronic  devices  has  recently  generated  considerable  interest  in  both  the  basic  transport 
physics  and  potential  technological  applications  of  “molecular  electronics”.3,4  However 
some  reports  of  molecular  mechanisms  in  electronic  devices5,6  have  been  shown  to  be 
premature  and  due  to  filamentary  conduction,7  highlighting  the  fabrication  sensitivity  of 
molecular  structures  and  the  need  to  institute  reliable  controls  and  methods  to  validate 
true  molecular  transport.8  A  related  problem  is  the  characterization  of  molecules  in  the 
active  device  structure,  including  their  configuration,  bonding,  and  indeed  even  their  very 
presence.  Here  we  present  results  on  well-understood  molecular  assemblies,  which 
exhibit  an  understood  classical  transport  behavior,  and  which  can  be  used  as  a  control  for 
eliminating  (or  understanding)  fabrication  variables.  Utilizing  tunneling  spectroscopic 
methods,  we  present  the  first  unambiguous  evidence  of  the  presence  of  molecules  in  the 
junction,  and  further  confirm  the  charge  transport  mechanism  obtained  by  standard 
current-voltage  characterizations. 

A  molecular  system  whose  structure  and  configuration  are  sufficiently  well- 
characterized  such  that  it  can  serve  as  a  standard  is  the  extensively  studied  alkanethiol 
(CH3(CH2)n-iSH)  self-assembled  monolayer  (SAM).9  This  system  is  useful  as  a  control 
since  properly  prepared  SAMs  form  single  van  der  Waals  crystals,9,10  and  presents  a 
simple  classical  metal-insulator-metal  (M-I-M)  tunnel  junction  when  fabricated  between 
metallic  contacts  due  to  the  large  HOMO-LUMO  gap  (HOMO:  highest  occupied 
molecular  orbital,  LUMO:  lowest  unoccupied  molecular  orbital)  of  approximately  8 
eV.11'13 

Various  surface  analytical  tools  have  been  utilized  to  investigate  the  surface  and  bulk 
properties  of  the  alkanethiol  SAMs,  such  as  X-ray  photoelectron  spectroscopy,14  Fourier 
transform  infrared  spectroscopy  (FTIR),15  Raman  spectroscopy,16  scanning  tunneling 
microscopy  (STM)10  etc.  For  example,  studies  have  shown  that  the  bonding  of  the 
thiolate  group  to  the  gold  surface  is  strong  with  a  bonding  energy  of  ~  1 .7  eV.9  STM 
topography  examinations  revealed  that  alkanethiols  adopt  the  commensurate  crystalline 
lattice  characterized  by  a  c(4  x  2)  superlattice  of  a  (V3  x  V3)R30°. 10,17  FTIR  investigation 
showed  that  the  orientation  of  the  alkanethiol  SAMs  on  Au(l  11)  surfaces  are  tilted  ~  30° 
from  the  surface  normal.18 

Electronic  transport  through  alkanethiol  SAMs  have  also  been  characterized  by 
STM, 19,20  conducting  atomic  force  microscopy,21'24  mercury-drop  junctions,25'28  cross- 

00  on 

wire  junctions,  and  electrochemical  methods.  '  These  investigations  are  exclusively 
at  ambient  temperature  -  clearly  useful  -  but  insufficient  for  an  unambiguous  claim  that 
the  transport  mechanism  is  tunneling  (of  course  expected,  assuming  that  the  Fermi  levels 
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of  the  contacts  lie  within  the  large  HOMO-LUMO  gap).  However  in  the  absence  of 
temperature-dependent  current- voltage  (I(V,T))  characteristics,  other  conduction 
mechanisms  (such  as  thermionic,  hopping,  or  filamentary  conduction)  cannot  be  excluded 
complicate  the  analysis,  and  thus  such  a  claim  is  premature. 

Utilizing  a  nanometer  scale  device  structure  that  incorporates  alkanethiol  SAMs,  we 
demonstrate  devices  that  allow  I(V,T)  and  structure-dependent  measurements33,34  with 
results  that  can  be  compared  with  accepted  theoretical  models  of  M-I-M  tunneling.  The 
use  of  this  fabrication  approach  is  not  special  in  any  way  (other  than  that  we  have  so  far 
found  it  to  be  successful)  -  indeed  we  stress  that  any  successful  device  fabrication 
method  should  yield  the  results  described  below  if  one  is  characterizing  the  intrinsic 
molecular  transport  properties. 

The  electronic  transport  is  further  investigated  with  the  technique  of  inelastic 
electron  tunneling  spectroscopy  (IETS).34  IETS  was  developed  in  the  1960’s  as  a 
powerful  spectroscopic  tool  to  study  the  vibrational  spectrum  of  organic  molecules 
confined  inside  metal-oxide-metal  junctions.35'39  In  our  study  IETS  is  utilized  for  the 
purpose  of  molecule  identification,  chemical  bonding  and  conduction  mechanism 
investigations  of  the  “control”  SAMs.  The  exclusive  presence  of  well-known  vibrational 
modes  of  the  alkanes  used  are  direct  evidence  of  the  molecules  in  the  device  structure, 
something  that  has  to  date  only  been  inferred  (with  good  reason,  but  nonetheless  not 
unambiguously).  The  vibrational  modes,  exclusively  identified  as  alkanes  (as  well  as 
contact  modes)  are  difficult  to  interpret  in  any  other  way  other  than  as  components  in  the 
active  region  of  the  device.  The  inelastic  tunneling  spectra  also  demonstrate  that 
electronic  tunneling  occurs  through  the  molecules,  confirming  the  conduction  mechanism 
obtained  by  I(V,T)  characterizations.  The  specific  spectral  lines  also  yield  intrinsic 
linewidths  that  may  give  insight  into  molecular  conformation,  and  may  prove  to  be  a 

powerful  tool  in  future  molecular  device 
characterization. 

Electronic  transport  measurements  on 
alkanethiol  SAMs  were  performed  using  a 
device  structure  similar  to  one  reported 
previously.33'34,40'42  In  this  device,  as  illustrated 
in  Figure  1(a)  (not  drawn  to  scale  in  the  relative 
thickness),  a  number  of  molecules  (~  several 
thousands)  are  sandwiched  between  two  metallic 
contacts.  This  technique  provides  a  stable  device 
structure  and  makes  cryogenic  measurements 
possible.  The  device  fabrication  starts  with  a 
high  resistivity  silicon  wafer  with  low  stress 
SbN4  film  deposited  on  both  sides  by  low 
pressure  chemical  vapor  deposition  (LPCVD). 

By  standard  photolithography  processing,  a 
suspended  SpN4  membrane  (size  of  40  pm  x  40 
pm  and  thickness  of  ~  70  nm)  is  fabricated  on 
the  topside  of  the  wafer.  Subsequent  e-beam 
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octanedithiol 


lithography  and  reactive  ion  etching  creates  a  single  pore  with  a  diameter  of  tens  of 
nanometers  through  the  membrane.  As  the  next  step,  150  nm  gold  is  thermally 
evaporated  onto  the  topside  of  the  wafer  to  fill  the  pore  and  form  one  of  the  metallic 
contacts. 


Figure  1.  Schematics  of  a  nanometer  scale  device  used  in  this  study  (not  drawn  to  scale  in  the 
relative  thickness),  (a)  Top  schematic  is  the  cross  section  of  a  silicon  wafer  with  a  nanometer 
scale  pore  etched  through  a  suspended  silicon  nitride  membrane.  Middle  and  bottom  schematics 
show  a  Au/SAM/Au  junction  formed  in  the  pore  area,  (b)  The  structures  of  octanethiol  and 
octanedithiol  are  shown  as  examples. 


The  device  is  then  transferred  into  a  molecular  solution  to  deposit  the  SAM  layer. 

For  our  experiments,  a  ~  5  mM  alkanethiol  solution  is  prepared  by  adding  ~  10  pL 
alkanethiols  into  10  mL  ethanol.43  The  deposition  is  done  in  solution  for  24  hours  inside  a 
nitrogen  filled  glove  box  with  an  oxygen  level  of  less  than  100  ppm.  Three 
alkanemonothiol  molecules  of  different  molecular  lengths-octanethiol  [CF^CF^SH; 
denoted  as  C8,  for  the  number  of  alkyl  units],  dodecanethiol  [CH3(CH2)nSH,  denoted  as 
C12],  and  hexadecanethiol  [CF^CEhjisSH,  denoted  as  C16]  and  one  alkanedithiol 
molecule-octanedithiol  [HS(CH2)gSH,  denoted  as  C8-dithiol]  were  used  to  form  the 
active  molecular  components.43  As  representative  examples,  the  chemical  structures  of 
octanethiol  and  octanedithiol  are  shown  in  Figure  1(b). 

In  order  to  statistically  determine  the  pore  size,  test  patterns  (arrays  of  pores)  were 
created  under  similar  fabrication  conditions.  Figure  2  shows  a  scanning  electron 
microscope  (SEM)  image  of  such  test  pattern  arrays.  This  indirect  measurement  of  device 
size  is  done  since  SEM  examination  of  the  actual  device  can  cause  hydrocarbon 
contamination  of  the  device  and  subsequent  contamination  of  the  monolayer.  From 
regression  analysis  of  298  pores,  the  device  sizes  of  the  C8,  Cl  2,  Cl  6,  and  C8-dithiol 

in  diameters, 
that  of  previous 
within  a  standard  error. 
We  will  use  these  device 
areas  as  the  effective 
contact  areas.  Although 
one  could  postulate  that 
the  actual  area  of  metal 
that  contacts  the 
molecules  may  be 
different,  there  is  little 
reason  to  propose  it 
would  be  different  as  a 
function  of  length  over 
the  range  of  alkanethiols 
used,  and  at  most  would 


samples  are  determined  as  50  ±  8,  45  ±  2,  45  ±  2,  and  5 1  ±  5  nm 
respectively.  A  more  ideal  (less  parasitic)  C8  sample  supercedes 
reports,33  and  derived  parameters  from  the  two  data  sets  agree  to 
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be  a  constant  systematic  error. 

Figure  2.  A  scanning  electron  microscope  image  of  a  representative  array  of  pores  used  to 
calibrate  device  size.  The  scale  bar  is  500  nm. 


The  sample  is  then  transferred  in  ambient  conditions  to  an  evaporator  that  has  a 
cooling  stage  to  deposit  the  opposing  Au  contact.  During  the  thermal  evaporation  (under 
the  pressure  of  ~  10'8  Torr),  liquid  nitrogen  is  kept  flowing  through  the  cooling  stage  in 
order  to  avoid  thermal  damage  to  the  molecular  layer.33,44  This  technique  reduces  the 
kinetic  energy  of  evaporated  Au  atoms  at  the  surface  of  the  monolayer,  thus  preventing 
Au  atoms  from  punching  through  the  monolayer.  For  the  same  reason  the  evaporation 
rate  is  kept  very  low.  For  the  first  10  nm  gold  evaporated,  the  rate  is  less  than  0.1  A/s. 
Then  the  rate  is  increased  slowly  to  0.5  A/s  for  the  rest  of  the  evaporation  and  a  total  of 
200  nm  gold  is  deposited  to  form  the  contact. 

The  device  is  subsequently  packaged  and  loaded  into  a  low  temperature  cryostat.  The 
sample  temperature  is  varied  from  300  to  4.2  K  by  flowing  cryogen  vapor  onto  the 
sample  (and  thermometer)  using  a  closed  loop  temperature  controller.  Two-terminal  dc 
I(V)  measurements  are  performed  using  a  semiconductor  parameter  analyzer.  Inelastic 
electron  tunneling  spectra  are  obtained  via  a  standard  lock-in  second  harmonic 
measurement  technique.35,36  A  synthesized  function  generator  is  used  to  provide  both  the 
modulation  and  the  lock-in  reference  signal.  The  second  harmonic  signal  (proportional  to 
d2I/dV2)  is  directly  measured  using  a  lock-in  amplifier,  which  is  checked  to  be  consistent 
with  a  numerical  derivative  of  the  first  harmonic  signal  (proportional  to  dl/dV).  Various 
modulation  amplitudes  and  frequencies  are  utilized  to  obtain  the  spectra.  The  ac 
modulation  is  added  to  a  dc  bias  using  operational  amplifier-based  custom  circuitry  45 

In  Table  1,  possible  conduction  mechanisms  are  listed  with  their  characteristic 
current,  temperature-  and  voltage-dependencies46  (We  do  not  discuss  filamentary 
tunneling  mechanisms,  which  are  easier  to  categorize47).  Based  on  whether  thermal 
activation  is  involved,  the  conduction  mechanisms  fall  into  two  distinct  categories:  (i) 
thermionic  or  hopping  conduction  which  has  temperature-dependent  I(V)  behavior  and 
(ii)  direct  tunneling  or  Fowler-Nordheim  tunneling  which  does  not  have  temperature- 
dependent  I(V)  behavior.  For  example,  thermionic  and  hopping  conductions  have  been 
observed  for  4-thioacetylbiphenyl  SAMs40  and  1, 4-phene lyene  diisocyanide  SAMs.4"3  On 
the  other  hand,  the  conduction  mechanism  is  expected  to  be  tunneling  when  the  Fermi 
levels  of  contacts  lie  within  the  large  HOMO-LUMO  gap  for  short  length  molecule,  as 
for  the  case  of  alkanethiol  molecular  system.11'13  Previous  work  on  Langmuir-Blodgett 
alkane  monolayers48  exhibited  a  significant  impurity-dominated  transport  component, 
complicating  the  analysis.  I(V)  measurements  on  self-assembled  alkanethiol  monolayers 
have  also  been  reported;19'29,49  however  all  of  these  measurements  were  performed  at 
fixed  temperature  (300  K)  which  is  insufficient  to  prove  tunneling  as  the  dominant 
mechanism. 
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TABLE  1.  Possible  conduction  mechanisms.  Adapted  from  Ref.  46. 


Conduction 

Mechanism 

Characteristic 

Behavior 

Temperature 

Dependence 

Voltage 

Dependence 

Direct 

Tunneling* 

J~Hexp(-— v/2wO) 
h 

none 

J~V 

Fowler- 

Nordheim 

Tunneling 

,  „2  ,  WST®3'2, 
J~Fexp(-  3  qnv  } 

none 

ln( — : r) - 

V 2  V 

Thermionic 

Emission 

kT 

1  1 

Wt2)~t 

1 

ln(J)  ~  V~2 

Hopping 

Conduction 

J  ~V  exp(  kT) 

ln(-)  ~  - 
V  T 

J~V 

*  This  characteristic  of  direct  tunneling  is  valid  for  the  low  bias  regime  [see  Eq.  (3a)]. 


To  describe  the  transport  through  a  molecular  system  having  HOMO  and  LUMO 
energy  levels,  one  of  the  applicable  models  is  the  Franz  two-band  model.50"53  This  model 
provides  a  non-parabolic  energy-momentum  E(k)  dispersion  relationship  by  considering 
the  contributions  of  both  the  HOMO  and  LUMO  energy  levels:50 


k2 


2jf *('♦ 


(1) 


where  k  is  the  imaginary  part  of  wave  vector  of  electrons,  m*  is  the  electron 
effective  mass,  h  (=  2nh)  is  Planck’s  constant,  E  is  the  electron  energy,  and  Eg  is  the 
HOMO-LUMO  energy  gap.  From  this  non-parabolic  E(k)  relationship,  the  effective  mass 
of  the  electron  tunneling  through  the  SAM  can  be  deduced  by  knowing  the  barrier  height 
of  the  metal-SAM-metal  junction. 
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When  the  Fermi  level  of  the  metal  is  aligned  close  enough  to  one  energy  level  (either 
HOMO  or  LUMO),  the  effect  of  the  other  distant  energy  level  on  the  tunneling  transport 
is  negligible,  and  the  widely  used  Simmons  model54  is  an  excellent  approximation.55 
Simmons  model  expressed  the  tunneling  current  density  through  a  barrier  in  the  tunneling 
regime  of  V  <  0B/e  as25,54 


J  = 


exp 


exp 


2(2m)' 


2(2  m)' 
h 


a  O , 


eV 

2  j 


(2) 


-a \  ®  B  + 


eV 


where  m  is  the  electron  mass,  d  is  the  barrier  width,  Ob  is  the  barrier  height,  and  V  is 
the  applied  bias.  For  molecular  systems,  the  Simmons  model  has  been  modified  with  a 
parameter  a.25,33  a  is  a  unitless  adjustable  parameter  that  is  introduced  to  provide  either  a 
way  of  applying  the  tunneling  model  of  a  rectangular  barrier  to  tunneling  through  a 
nonrectangular  barrier,25  or  an  adjustment  to  account  for  the  effective  mass  (m*)  of  the 
tunneling  electrons  through  a  rectangular  barrier,  ’  ’  ’  or  both,  a  =  1  corresponds  to 
the  case  for  a  rectangular  barrier  and  bare  electron  mass.  By  fitting  individual  I(V)  data 
using  Eq.  (2),  Ob  and  a  values  can  be  obtained. 

Eq.  (2)  can  be  approximated  in  two  limits:  low  bias  and  high  bias  as  compared  with 
the  barrier  height  Ob.  For  the  low  bias  range,  Eq.  (2)  can  be  approximated  as54 
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(2m<&  B)'ne2a 


h  d 


V  exp 


n 


(3a) 


To  determine  the  high  bias  limit,  we  compare  the  relative  magnitudes  of  the  first  and 
second  exponential  terms  in  Eq.  (2).  At  high  bias,  the  first  term  is  dominant  and  thus  the 
current  density  can  be  approximated  as 
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(3b) 


The  tunneling  currents  in  both  bias  regimes  are  exponentially  dependent  on  the 

barrier  width  d.  In  the  low  bias  regime  the  tunneling  current  density  is  J  oc  —  exp (~(3ad) , 

d 

where  Po  is  bias-independent  decay  coefficient: 


Po  = 


2(2  m)U2 
h 


(4a) 
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while  in  the  high  bias  regime,  J  oc  —exp (-/3vd) ,  where  (3V  is  bias-dependent  decay 


coefficient: 


effect  due  to  the  applied  bias. 

In  order  to  determine  the  conduction  mechanism  of  self-assembled  alkanethiol 
molecular  systems  I(V)  measurements  in  a  sufficiently  wide  temperature  range  (300  to  80 
K)  and  resolution  (10  K)  were  performed.  Figure  3(a)  shows  a  representative  I(V,T) 
characteristic  of  dodecanethiol  (Cl 2)  measured  with  the  device  structure  as  shown  in 
Figure  1(a).  Positive  bias  corresponds  to  electrons  injected  from  the  physisorbed  Au 
contact  [bottom  contact  in  Figure  1(a)]  into  the  molecules.  By  using  the  contact  area  of 
45  ±  2  nm  in  diameter  determined  from  SEM  study,  a  current  density  of  1,500  ±  200 
A/cm2  at  1.0  Volt  is  determined.  No  significant  temperature  dependence  of  the 
characteristics  (from  V  =  0  to  1.0  Volt)  is  observed  over  the  range  from  300  to  80  K.  An 
Arrhenius  plot  (ln(I)  versus  1/T)  of  this  is  shown  in  Figure  3(b),  exhibiting  little 
temperature  dependence  in  the  slopes  of  ln(I)  versus  1/T  at  different  bias  and  thus 
indicating  the  absence  of  thermal  activation.  Therefore,  we  conclude  that  the  conduction 
mechanism  through  alkanethiol  is  tunneling  contingent  on  demonstrating  a  correct 
molecular  length  dependence.  The  tunneling  through  alkanethiol  SAMs  has  been 
assumed  as  “through-bond”  tunneling,  i.e.,  along  the  tilted  molecular  chains  between  the 
metal  contacts.  ’  ’  ’  Based  on  the  applied  bias  as  compared  with  the  barrier  height 
(Ob),  the  tunneling  through  a  SAM  layer  can  be  categorized  into  either  direct  (V  <  0B/e) 
or  Fowler-Nordheim  (V  >  Oe/e)  tunneling.  These  two  tunneling  mechanisms  can  be 
distinguished  due  to  their  distinct  voltage  dependencies  (see  Table  1).  Analysis  of 
ln(I/V2)  versus  1/V  [in  Figure  3(c)]  shows  no  significant  voltage  dependence,  indicating 
no  obvious  Fowler-Nordheim  transport  behavior  in  this  bias  range  (0  to  1.0  Volt)  and 
thus  determining  that  the  barrier  height  is  larger  than  the  applied  bias,  i.e.,  Ob  >  1.0  eV. 
This  study  is  restricted  to  applied  biases  <  1 .0  Volt  and  the  transition  from  direct  to 
Fowler-Nordheim  tunneling  requires  higher  bias. 
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Figure  3.  (a) 
Temperature- 
dependent  I(V) 
characteristics  of 
dodecanethiol  (Cl 2). 
I(V)  data  at 
temperatures  from 
300  to  80  K  with  20 
K  steps  are  plotted 
on  a  log  scale,  (b) 
Arrhenius  plot 
generated  from  the 
I(V)  data  in  (a),  at 
voltages  from  0.1  to 
1.0  Volt  with  0.1 
Volt  steps,  (c)  Plot  of 
ln(I/V2)  versus  1/V 
at  selected 
temperatures. 


The 

importance  of 
variable 
temperature 
measurements  to 
validate  tunneling 
is  demonstrated  in 
Figure  4.  Here  the 
I(V)  of  an 
octanethiol  (C8) 
device  is  shown 
(Figure  4(a)), 
whose  I(V)  shape 
looks  very  similar 
to  Figure  3  (i.e., 
direct  tunneling), 
and  indeed  can  be 
fit  to  a  Simmons 
model.  However, 
further  I(V,T) 
measurements 
display  an  obvious 
temperature 
dependence 
(Figure  4(b)), 
which  can  be  fit 
well  to  a  hopping 
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conduction  model  (Table  1)  with  a  well-defined  activation  energy  of  190  meV  (Figure 

4(c)).  This  and  other 
similar  impurity- 
mediated  transport 
phenomena  (such  as 
Coulomb  blockade)  are 
observed  in  a  subset  of 
devices  and  is 
indicative  of  the 
unintentional 
incorporation  of  a  trap 
or  defect  level  in  those 
devices.  This  study 
instead  focuses  on 
devices  which  do  not 
show  any  defect- 
mediated  transport  and 
probes  the  intrinsic 
behavior  of  the 
molecular  layer. 


Figure  4.  (a)  I(V) 
characteristics  of  an 
octanethiol  (C8)  device  at 
270K.  (b)  Temperature 
dependence  of  the  device 
from  270K  to  1 80K  (in 
10K  increments),  (c)  Plot 
of  ln(I/V)  versus  1/T  at 
various  voltages.  The 
activated  behavior  is  bias 
voltage  independent,  thus 
the  behavior  is  hopping 
(in  this  device)  due  to 
incorporation  of  a  defect 
of  energy  1 90  meV.  This 
class  of  device  is  not 
suitable  for  investigation 
of  the  intrinsic  transport 
mechanism  in  the  SAM 
as  it  is  dominated  by  a 
defect. 

Having  established 
tunneling  as  the 

3.5  4.0  4.5  5.g8  5.5 

1000/T  (1/K) 


conduction  mechanism  in  a  device,  we  will  now  obtain  the  barrier  height  by  comparing 
experimental  I(V)  data  with  theoretical  calculations  from  tunneling  models. 

From  the  modified  Simmons  model  [Eq.  (2)]  by  adjusting  two  parameters  <DB  and  a, 
a  nonlinear  least  squares  fitting  can  be  performed  to  fit  the  measured  C12  I(V)  data 
(calculation  assuming  a  =  1  has  been  previously  shown  not  to  fit  I(V)  data  well  for  some 
alkanethiol  measurements  at  fixed  temperature  (300  K)).25  By  using  a  device  size  of  45 
nm  in  diameter,  the  best  fitting  parameters  (minimizing  y2)  for  the  room  temperature  C12 
I(V)  data  were  found  to  be  Ob  =  1.42  ±  0.04  eV  and  a  =  0.65  ±0.01,  where  the  error 
ranges  of  Ob  and  a  are  dominated  by  potential  device  size  fluctuations  of  2  nm. 

Likewise,  data  sets  were  obtained  and  fittings  were  done  for  octanethiol  (C8)  and 
hexadecanethiol  (C16),  which  yielded  values  {Ob  =  1.83  ±  0.10  eV  and  a  =  0.61  ±  0.01} 
and  {Ob  =  1.40  ±  0.03  eV,  a  =  0.68  ±  0.01},  respectively. 

Using  Ob  =  1.42  eV  and  a  =  0.65,  a  calculated  I(V)  for  C12  is  plotted  as  a  solid 
curve  in  a  linear  scale  [Figure  5(a)]  and  in  a  semi-log  scale  [Figure  5(b)].  A  calculated 
I(V)  for  a  =  1  and  Ob  =  0.65  eV  (which  gives  the  best  fit  at  low  bias  range)  is  shown  as 
the  dashed  curve  in  the  same  figure,  illustrating  that  with  a  =  1  only  limited  regions  of 
the  I(V)  can  be  fit  (specifically  here,  for  V  <  0.3  Volt).  For  the  case  of  a  rectangular 
barrier,  the  a  parameter  fit  presented  above  corresponds  to  an  effective  mass  m*  ( =  a2 


m)  of  0.42  m. 


Figure  5.  Measured  C12  I(V)  data  (circular  symbols)  is  compared  with  calculation  (solid  curve) 
using  the  optimum  fitting  parameters  of  <hB  =  1-42  eV  and  a  =  0.65.  The  calculated  I(V)  from  a 
simple  rectangular  model  (a  =  1)  with  4>B  =  0.65  eV  is  also  shown  as  the  dashed  curve.  Current  is 
plotted  (a)  on  linear  scale  and  (b)  on  log  scale. 
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In  order  to  investigate  the  dependency  of  the  Simmons  model  fitting  on  Ob  and  a,  a 
fitting  minimization  analysis  was  undertaken  on  the  individual  Ob  and  a  values  as  well  as 
their  product  form  of  ad>B1/2  in  Eq.  (4a).  A(Ob,  a)  =  (  E  |  Iexp>v  -  ICai,v  |2  )1/2  was  calculated 
and  plotted  where  IeXp,v  is  the  experimental  current-voltage  values  and  Icai,v  is  calculated 
using  Eq.  (2).  7,500  different  {Ob,  a}  pairs  were  used  in  the  fittings  with  Ob  ranging 
from  1.0  to  2.5  eV  (0.01  eV  increment)  and  a  from  0.5  to  1.0  (0.01  increment).  Figure 
6(a)  is  a  representative  contour  plot  of  A(Ob,  a)  versus  Ob  and  a  values  generated  for  the 
C12  I(V)  data  where  darker  regions  correspond  to  smaller  A(Ob,  a)  and  various  shades 
represent  half  order  of  magnitude  A(Ob,  a)  steps.  The  darker  regions  represent  better  fits 
of  Eq.  (2)  to  the  measured  I(V)  data.  In  the  inset  in  Figure  6(a)  one  can  see  there  is  a 
range  of  possible  Ob  and  a  values  yielding  minimum  fitting  parameters.  Although  the 
tunneling  parameters  determined  from  the  previous  Simmons  tunneling  fitting  {®B  = 

1.42  eV  and  a  =  0.65}  lie  within  this  minimum  region  in  this  figure,  there  is  a 
distribution  of  other  possible  values. 


03 
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Figure  6.  (a)  Contour  plot  of 
A(Ob,  a)  values  for  Cl 2 
nanopore  device  as  a  function  of 
Ob  and  a,  where  the  darker 
region  corresponds  to  a  better 
fitting.  Inset  shows  detailed 
minimization  fitting  regions,  (b) 
A  plot  of  A(Ob,  a)  as  a  function 
of  a  Ob1/2. 

A  plot  of  A(Ob,  a) 
versus  ad>B1/2  for  the  same 
device  reveals  a  more 
pronounced  dependence,  and 
is  shown  in  Figure  6(b).  This 
plot  indicates  the  fitting  to 
the  Simmons  model  sharply 
depends  on  the  product  of 
a<E>e1/2.  For  this  plot  the 
A(<E>b,  a.)  is  minimized  at 
aC>B1/2  of  0.77  (eV)1/2 
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corresponding  to  a  (30  value 
of  0.79  A"!  from  Eq.  (4a).  The  C8  and  C16 
devices  showed  similar  results,  indicating 
the  Simmons  tunneling  model  has  a  strong 
a<J>B  dependence.  For  the  C8  device, 
although  <1>b  obtained  from  the  fitting  is  a 


little  larger,  combined  a  and  <DB  gives  a  similar  p0  value  within  the  error  range  as  the  C12 
and  C16  devices  (Table  2). 


TABLE  2.  Summary  of  alkanethiol  tunneling  parameters  in  this  study 


Molecules 

J  at  1  V  (A/cm2) 

(eV) 

a 

m*  (m) 

Po(A-‘)* 

C8 

31,000  ±  10,000 

1.83  ±0.10 

0.61  ±0.01 

0.37 

0.85  ±  0.04 

C12 

1,500  ±200 

1.42  ±0.04 

0.65  ±0.01 

0.42 

0.79  ±  0.02 

C16 

23+2 

1.40  ±0.03 

0.68  ±0.01 

0.46 

0.82  ±  0.02 

C8-dithiol 

93,000  ±  18,000 

1.20  ±0.03 

0.59  ±0.01 

0.35 

0.66  ±  0.02 

*  Po  values  were  calculated  using  Eq.  (4a). 

Three  alkanethiols  of  different  molecular  length,  C8,  Cl  2,  and  C16  were  investigated 
to  study  length-dependent  tunneling  behavior.  Figure  7  is  a  semi-log  plot  of  tunneling 
current  densities  multiplied  by  molecular  length  (Jd  at  low  bias  and  Jd2  at  high  bias)  as  a 
function  of  the  molecular  length  for  these  alkanethiols.  The  molecular  lengths  used  in  this 
plot  are  13.3,  18.2,  and  23.2  A  for  C8,  C12,  and  C16,  respectively.  Each  molecular  length 
was  determined  by  adding  an  Au-thiol  bonding  length  to  the  length  of  molecule22.  Note 
that  these  lengths  assume  through-bond  tunneling.  ’  ’  '  The  high  and  low  bias  regimes 
are  defined  somewhat  arbitrarily  by  comparing  the  relative  magnitudes  of  the  first  and 
second  exponential  terms  in  Eq.  (2).  Using  Ob  =  1 .42  eV  and  a  =  0.65  obtained  from 
nonlinear  least  squares  fitting  of  the  C 12  I(V)  data,  the  second  term  becomes  less  than  ~ 
10  %  of  the  first  term  at  ~  0.5  Volt  that  is  chosen  as  the  boundary  of  low  and  high  bias 
ranges. 
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Figure  7.  Log  plot  of  tunneling  current  densities  multiplied  by  molecular  length  d  at  low  bias  and 
by  d2  at  high  bias  (symbols)  versus  molecular  lengths.  The  lines  through  the  data  points  are  linear 
fittings. 

As  seen  in  Figure  7,  the  tunneling  current  shows  exponential  dependence  on 
molecular  length,  which  is  consistent  with  the  Simmons  tunneling  model  [Eq.  (3)].  The  p 
values  can  be  determined  from  the  slope  at  each  bias  and  are  plotted  in  Figure  8.  The 
error  bar  of  an  individual  p  value  in  this  plot  was  obtained  by  considering  both  the  device 
size  uncertainties  and  the  linear  fitting  errors. 
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Jd2  (A) 


Figure  8.  Plot  of  p 
versus  bias  in  the 
low  bias  range 
(square  symbols) 
and  high  bias 
ranges  (circular 
symbols).  The  inset 
shows  a  plot  of  p2v 
versus  bias  with  a 
linear  fitting. 


The 

determined  p 
values  are  almost 
independent  of 
bias  in  the  low 


bias  range  (V  <  ~  0.5  V),  and  an  average  p  of  0.77  ±  0.06  A'1  in  this  region  (from  0  to  0.5 
V)  can  be  calculated  from  Figure  8.  Table  3  is  a  summary  of  previously  reported 
alkanethiol  transport  parameters  obtained  by  different  techniques.  The  current  densities 
(J)  listed  in  Table  3  are  for  C12  monothiol  or  dithiol  devices  at  1  V,  which  are 
extrapolated  from  published  results  of  other  length  alkane  molecules.  The  large  variation 
of  J  of  these  reports  can  be  attributed  to  the  uncertainties  in  device  contact  geometry  and 
junction  area,  as  well  as  complicating  inelastic  or  defect  contributions.  The  P  value  (0.77 
±  0.06  A"1  «  0.96  ±  0.08  per  methylene)  for  alkanethiols  reported  here  is  comparable  to 
previously  reported  values  as  summarized  in  Table  3.  This  p  value  agrees  with  the  value 
of  0.79  A’1  (Po)  calculated  via  Eq.  (4a)  from  fitting  individual  I(V)  characteristic  of  the 


C12  device.  The  calculated  p0  of  C8  and  Cl 6  devices  also  have  similar  values,  as 
summarized  in  Table  2. 


TABLE  3.  Summary  of  alkanethiol  tunneling  characteristic  parameters 


Junction 

P  (A'1) 

J  (A/cm2)  at  1  V 

G>b  (eV) 

Technique 

Ref. 

(bilayer)  monothiol 

0.87±0.1 

25-200a) 

2.1e) 

Hg-junction 

25 

(bilayer)  monothiol 

0.71±0.08 

0.7-3. 5a) 

Hg-junction 

27 

monothiol 

0.79±0,01 

1500±200b) 

1.4e) 

Solid  M-I-M 

33 

monothiol 

1.2 

STM 

19 

dithiol 

0.8+0.08 

3.7-5  x  105c) 

5  +  2° 

STM 

20 

monothiol 

0.73-0.95 

1 100-1900d) 

2.2e) 

CAFM 

21 
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monothiol 

0.64-0.8 

10-50d) 

2.3e) 

CAFM 

23 

dithiol 

0.46+0.02 

3-6  x  105c) 

1.3-1.5e) 

CAFM 

24 

monothiol 

1.37+0.03 

1.815 

Tuning  fork  AFM 

49 

monothiol 

0.97+0.04 

Electrochemical 

30 

monothiol 

0.85 

Electrochemical 

31 

monothiol 

0.91+0.08 

Electrochemical 

32 

monothiol 

0.76 

2xl04  (at  0.1  V)c) 

1.3-3.4®5 

Theory 

58 

monothiol 

0.76 

Theory 

59 

monothiol 

0.79 

Theory 

56 

Note: 

Some  decay  coefficients  p  were  converted  into  the  unit  of  A'1  from  the  unit  of  per  methylene. 

The  junction  areas  were  estimated  by  optical  microscope*5,  SEMb),  assuming  single  molecule65,  and 
Hertzian  contact  theoryd). 

Current  densities  (J)  for  C12  monothiol  or  dithiol  at  1  V  are  extrapolated  from  published  results  for  other 
length  molecules  by  using  conductance  °c  exp(-P  d)  relationship. 

Barrier  height  ®B  values  were  obtained  from  Simmons  equation65,  bias-dependence  of  fSf),  and  a  theoretical 
calculation85. 

According  to  Eq.  (4b),  p2v  depends  on  bias  V  linearly  in  the  high  bias  range.  The 
inset  in  Figure  8  is  a  plot  of  p2v  versus  V  in  this  range  (0.5  to  1 .0  Volt)  along  with  linear 
fitting  of  the  data.  From  this  fitting,  Ob  =  1 .35  ±  0.20  eV  and  a  =  0.66  ±  0.04  were 
obtained  from  the  intercept  and  the  slope,  respectively,  consistent  with  the  values  {Ob  = 
1.42  eV  and  a  =  0.65}  obtained  from  the  nonlinear  least  squares  fitting  in  the  previous 
section. 

P  values  for  alkanethiols  obtained  by  various  experimental  techniques  have 
previously  been  reported  and  are  summarized  in  Table  3. 19-33,49  In  order  to  compare  with 
these  reported  p  values,  we  also  performed  length-dependent  analysis  on  our 
experimental  data  according  to  the  generally  used  equation:20-28,33 

G  =  G0  exp(-/?c/)  (5) 

This  gives  a  p  value  from  0.84  to  0.73  A-1  in  the  bias  range  from  0.1  to  1.0  volt, 
which  is  comparable  to  results  reported  previously.  For  example,  Holmlin,  et.  al,  reported 
a  p  value  of  0.87  A-1  by  mercury  drop  experiments,25  and  Wold,  et.  al,  have  reported  P  of 
0.94  A-1  and  Cui,  et.  al,  reported  P  of  0.64  A-1  for  various  alkanethiols  by  using  a 
conducting  atomic  force  microscope  technique.21,23  These  reported  p  were  treated  as  bias- 
independent  quantities,  contrary  to  the  results  reported  here  and  that  observed  in  a 
slightly  different  alkane  system  (ligand-encapsulated  nanoparticle/alkane-dithiol 
molecules).24  We  also  caution  against  the  use  of  parameters  that  have  not  been  checked 
with  a  temperature-dependent  analysis,  since  small  non-tunneling  components  can 
dramatically  affect  derived  values  of  p. 
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We  have  analyzed  our  experimental  data  using  a  Franz  two-band  model.50'53  Since 
there  is  no  reliable  experimental  data  on  the  Fermi  level  alignment  in  these  metal-SAM- 
metal  systems,  <E>b  and  m*  are  treated  as  adjustable  parameters.  We  performed  a  least 
squares  fit  on  our  data  with  the  Franz  non-parabolic  E(k)  relationship  [Eq.  (1)]  using  an 
alkanethiol  HOMO-LUMO  gap  of  8  eV. 12,13  Figure  9  shows  the  resultant  E(k) 
relationship  and  the  corresponding  energy  band  diagrams.  The  zero  of  energy  in  this  plot 
was  chosen  as  the  LUMO  energy.  The  best  fitting  parameters  obtained  by  minimizing  %2 
were  Ob  =  1 .49  ±  0.5 1  eV  and  m*  =  0.43  ±  0. 1 5  m,  where  the  error  ranges  of  Ob  and  m* 
are  dominated  by  the  error  fluctuations  of  p  [k2  =  -  (p/2)2].  Both  electron  tunneling  near 
the  LUMO  and  hole  tunneling  near  the  ElOMO  can  be  described  by  these  parameters.  Ob 

=  1 .49  eV  indicates  that 
the  Fermi  level  is  aligned 
close  to  one  energy  level 
in  either  case,  therefore 
the  Simmons  model  is  a 
valid  approximation.  The 
Ob  and  m*  values 
obtained  here  are  in 
reasonable  agreement  with 
the  previous  results 
obtained  from  the 
Simmons  model. 


Figure  9.  E(k)  relationship 
(symbols)  generated  from  the 
length-dependent 
measurement  data  for 
alkanethiols.  Solid  and  open 
symbols  correspond  to 
electron  and  hole  tunneling, 
respectively.  The  insets  show 
the  corresponding  energy 
band  diagrams.  The  solid 
curve  is  the  Franz  two-band 


Electronic  transport  through  alkanethiol  SAMs  is  further  investigated  with  the 
technique  of  inelastic  electron  tunneling  spectroscopy,34  such  as  the  works  of  1966  by 
Jaklevic  and  Lambe  who  studied  the  conductance  of  a  tunnel  junctions  with  encased 
organic  molecules.35  Since  then  it  has  become  a  powerful  spectroscopic  tool  for  chemical 
identification,  chemical  bonding  investigation,  and  surface  chemistry  and  physics 
studies.38  In  an  inelastic  tunneling  process  the  electron  loses  energy  to  a  localized 
vibrational  mode  with  a  frequency  v  when  the  applied  bias  satisfies  the  condition  of  eV  = 
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hv.  As  a  result,  an  additional  tunneling  channel  is  opened  for  the  electron,  resulting  in  an 
increase  in  the  total  current  at  the  applied  bias  corresponding  to  the  vibrational  mode 
energy.37  Typically  only  a  small  fraction  of  tunneling  electrons  are  involved  in  the 
inelastic  tunneling  process  (determined  by  the  electron  -  vibronic  mode  coupling 
coefficient),  resulting  in  a  small  conductance  change,  which  is  commonly  measured  in 
the  second  harmonics  of  a  phase-sensitive  detector  that  yields  the  characteristic 
frequencies  of  the  corresponding  vibrational  modes  as  well  as  other  information.36'38 


I(V,T)  measurements  and  additional  IETS  studies  have  been  performed  on  an 
octanedithiol  (C8-dithiol)  SAM  using  the 
aforementioned  device  structure  shown  in 
Figure  1(a).34  Figure  10(a)  is  the  I(V,T)  data 
for  this  device  obtained  from  300  to  4.2  K. 

An  Arrhenius  plot  shown  in  Figure  10(b) 
exhibits  little  temperature  dependence, 
verifying  that  tunneling  is  the  main  transport 
mechanism  for  C8-dithiol  SAM.  This  result  is 
in  good  agreement  with  the  tunneling 
transport  characteristics  observed  previously. 

Figure  10(c)  shows  the  room  temperature  I(V) 
measurement  result.  Using  a  junction  area  of 
5 1  ±  5  nm  in  diameter  (obtained  from 
statistical  studies  of  the  nanopore  size  with 
SEM),  a  current  density  of  (9.3  ±  1.8)  x  104 
A/cm2  at  1.0  Volt  is  calculated.  As  a 
comparison,  the  current  density  of  (3.1  ±  1.0) 
x  104  A/cm2  at  1.0  Volt  was  observed  for  C8 
monothiol  SAM.  Using  the  modified 
Simmons  model  [Eq.  (2)],  the  transport 
parameters  of  <J>B  =  1 .20  ±  0.03  eV  and  a  = 

0.59  ±  0.01  (m*  =  0.34  m)  were  obtained  for 
this  C8-dithiol  SAM. 


Figure  10.  (a)  I(V,T)  characteristics  of  C8  dithiol 
SAM  at  selected  temperatures  (4.2,  50,  100,  150, 
200,  250,  and  290  K).  (b)  Arrhenius  plot 
generated  from  the  data  in  (a),  at  voltages  from 
0.1  to  0.5  Volt  with  0.05  Volt  steps,  (c)  Measured 
C8-dithiol  I(V)  data  at  room  temperature  (circular 
symbols)  is  compared  with  calculation  (solid 
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curve)  using  the  optimum  fitting  parameters  of  <&B  =  1 .20  eV  and  a  =  0.59. 


Figure  1 1  shows  the  IETS  spectrum  of  the  same  C8-dithiol  SAM  device  obtained  at 
T  =  4.2  K.  An  ac  modulation  of  8.7  mV  (rms  value)  at  a  frequency  of  503  Hz  was  applied 
to  the  sample  to  acquire  the  second  harmonic  signals.  The  spectra  are  stable  and 
repeatable  upon  successive  bias  sweeps.  The  spectrum  at  4.2  K  is  characterized  by  three 
pronounced  peaks  in  the  0  to  200  mV  region  at  33,  133,  and  158  mV.  From  comparison 
with  previously  reported  infrared  (IR),  Raman,  and  high  resolution  electron  energy  loss 
(HREEL)  spectra  of  SAM  covered  gold  surfaces  (Table  4),  these  three  peaks  are  assigned 
to  v(Au-S),  v(C-C),  and  yw(CH2)  modes  of  a  surface  bound  alkanethiolate.60'63  The 
absence  of  a  strong  v(S-H)  signal  at  -329  mV  suggests  that  most  of  the  thiol  groups  have 
reacted  with  the  gold  bottom  and  top  contacts.  Peaks  are  also  reproducibly  observed  at 
80,  107  mV,  and  186  mV.  They  correspond  to  v(C-S),  5r(CH2),  and  8S(CH2)  modes.  The 
stretching  mode  of  the  CH2  groups,  v(CH2),  appears  as  a  shoulder  at  357  meV.  The  peak 
at  15  mV  is  due  to  vibrations  from  either  Si,  Au,  or  5(C-C-C).64  We  note  that  all 
alkanethiolate  peaks  without  exception  or  omission  occur  in  the  spectra.  Peaks  at  58,  257, 
277,  and  302,  as  well  as  above  375  mV  are  likely  to  originate  from  Si-H  and  N-H 
vibrations  related  to  the  silicon  nitride  membrane, 643,65  which  forms  the  SAM 
encasement.  To  the  best  of  our  knowledge  aikanethiols  have  no  vibrational  signatures  in 
these  regions.  Measurement  of  the  background  spectrum  of  an  “empty”  nanopore  device 
with  only  gold  contacts  to  obtain  background  contributions  from  Si3N4  is  hampered  by 
either  too  low  (open  circuit)  or  too  high  (short  circuit)  currents  in  such  a  device.  Similar 
IETS  result  has  also  been  obtained  using  a  different  test  structure  recently.66 
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Figure  11.  Inelastic  electron  tunneling  spectrum  of  C8  dithiol  SAM  obtained  from  lock-in  second 
harmonic  measurement  with  an  ac  modulation  of  8.7  mV  (rms  value)  at  a  frequency  of  503  Hz  (T 
=  4.2  K).  Peaks  labeled  with  Si  are  most  probably  background  due  to  the  encasing  S^N^ 


TABLE  4:  Summary*  of  the  major  vibrational  modes  of  alkanethiolates.  Taken  from  Ref. 
60-62. 


Modes 

Methods 

Wavenumber  (cm'1) 

(meV) 

v(Au-S) 

HREELS61 

225 

28 

Raman60 

641 

79 

Raman60 

706 

88 

HREEL61 

715 

89 

IR62 

720 

89 

IR62 

766 

95 

IR62 

925 

115 

HREEL61 

1050 

130 

78 


v(C-C) 

Raman60 

1064 

132 

Raman60 

1120 

139 

IR62 

1230 

152 

HREELS61 

1265 

157 

IR62 

1283 

159 

IR62 

1330 

165 

5s(CH2) 

HREELS61 

1455 

180 

v(S-H) 

Raman60 

2575 

319 

Raman60 

2854 

354 

HREELS61 

2860 

355 

Raman60 

2880 

357 

Raman60 

2907 

360 

HREELS61 

2925 

363 

*  There  is  a  vast  amount  of  literature  with  spectroscopic  assignments  for  alkanethiols. 
The  references  given  are  representative  for  IR,  Raman,  and  HREELS  assignments. 


Although  there  are  no  selection  rules  in  IETS  as  there  are  in  IR  and  Raman 
spectroscopy,  certain  selection  preferences  have  been  established.  According  to  the  IETS 
theory,67  molecular  vibrations  with  net  dipole  moments  perpendicular  to  the  interface  of 
the  tunneling  junction  have  stronger  peak  intensities  than  vibrations  with  net  dipole 
moments  parallel  to  the  interface  (for  dipoles  close  to  the  electrodes).  Thus  vibrations 
perpendicular  to  the  electrode  interface,  i.e.,  v(Au-S),  v(C-S),  v(C-C),  and  yw(CH2), 
dominate  the  IETS  spectrum  while  modes  parallel  to  the  interface,  i.e.,  5r,s(CH2)  and 
v(CH2),  are  weak,  as  clearly  shown  in  Figure  1 1 . 

In  order  to  verify  that  the  observed  spectra  are  indeed  valid  IETS  data,  the  peak 
width  broadening  was  examined  as  a  function  of  temperature  and  modulation  voltage. 
IETS  was  performed  with  different  ac  modulations  at  a  fixed  temperature,  and  at  different 
temperatures  with  a  fixed  ac  modulation.  Figure  12(a)  shows  the  modulation  dependence 
of  the  IETS  spectra  obtained  at  4.2  K,  and  Figure  12(b)  shows  the  modulation  broadening 
of  the  C-C  stretching  mode  at  133  meV.  The  circular  symbols  are  the  full  widths  at  half 
maximum  (FWHMs)  of  the  experimental  peak  at  T  =  4.2  K  with  various  modulation 
voltages.  A  Gaussian  distribution  function  was  utilized  to  obtain  a  FWHM  and  the  error 
range.68  The  square  symbols  are  calculated  FWHM  values  (Wtheoreticai)  taking  into  account 
both  a  finite  temperature  effect  (Wthermai  ~  5.4  kBT)36  and  a  finite  voltage  modulation 
effect  (Wmoduiation  ~  1.7  Vac_rms)-69  These  two  broadening  contributions  add  as  the  squares: 
W2theoreticai  =  W2thcrmai  +  W2moduiation-  The  agreement  is  excellent  over  most  of  the 
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modulation  range,  but  we  note  a  saturation  of  the  linewidth  at  low  modulation  bias 
indicating  the  influence  of  a  non-negligible  intrinsic  linewidth.  Taking  into  account  the 
known  thermal  and  modulation  broadenings,  and  including  the  intrinsic  linewidth  (Wi)70 
as  a  fitting  parameter,  the  measured  peak  width  (Wexp)  is  given  by 

(6) 

Wi  can  be  determined  by  using  a  nonlinear  least  squares  fit  to  the  ac  modulation  data 
(Figure  12)  with  Eq.  (6),  giving  an  intrinsic  linewidth  of  3.73  ±  0.98  meV  for  this  line. 
This  is  shown  (with  the  error  range)  in  Figure  12(b)  as  a  shaded  bar,  including  the 
thermal  contribution. 
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Figure  12.  (a)  Modulation  dependence  of 
IETS  spectra  obtained  at  4.2  K.  (b)  Line 
(C-C  stretching  mode)  broadening  as  a 
function  of  ac  modulation.  The  circular 
symbols  are  experimental  FWHMs  and 
the  square  symbols  are  theoretical 

calculations  considering  both  modulation  and  thermal  contributions.  The  shaded  bar  denotes  the 
expected  saturation  due  to  the  derived  intrinsic  linewidth  (including  a  5.4  kBT  thermal 
contribution)  of  3.73  ±  0.98  meV. 


We  can  independently  check  the  thermal  broadening  of  the  line  at  fixed  modulation. 
Figure  13(a)  shows  the  temperature  dependence  of  the  IETS  spectra  obtained  with  an  ac 
modulation  of  8.7  mV  (rms  value).  In  Figure  13(b)  the  circular  symbols  (and 
corresponding  error  bars)  are  experimental  FWHM  values  of  the  C-C  stretching  mode 
from  Figure  13(a),  determined  by  a  Gaussian  fit  (and  error  of  the  fit)  to  the  experimental 
lineshape.  For  simplicity  we  have  only  considered  Gaussian  lineshapes68  resulting  in 
increased  error  bars  for  the  lower  temperature  range  due  to  an  asymmetric  lineshape.  The 
square  symbols  are  theoretical  calculations  considering  thermal  broadening,  modulation 
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broadening,  and  the  intrinsic  Iinewidth  determined  above.  The  error  ranges  of  the 
calculation  (due  to  the  intrinsic  Iinewidth  error)  are  approximately  the  size  of  the  data 
points.  The  agreement  between  theory  and  experiment  is  very  good,  spanning  a 
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temperature  range  from  below  (x  0.5) 
to  above  (xlO)  the  thermally 
broadened  intrinsic  Iinewidth.  This 
Iinewidth  should  be  a  sensitive  test  to 
compare  to  theoretical  models  of 
transmission  probabilities.71 


Figure  13.  (a)  Temperature  dependence  of 
IETS  spectra  obtained  at  a  fixed  ac  modulation  of  8.7  mV  (rms  value),  (b)  Line  (C-C  stretching 
mode)  broadening  as  a  function  of  temperature.  The  circular  symbols  are  experimental  FWHMs 
and  the  square  symbols  are  theoretical  calculations  considering  thermal  broadening,  modulation 
broadening,  and  the  intrinsic  Iinewidth.  . 


AC  modulation  (RMS  value)  (mV) 


Similar  intrinsic  linewidths  have  been 
determined  for  the  Au-S  stretching  mode  (33 
meV)  and  the  CH2  wagging  mode  (158  meV), 
and  the  results  are  shown  in  Figure.  14.  For 
the  Au-S  stretching  mode,  the  deviation  of 
experimental  data  from  calculated  values 
(thermal  and  modulation  width  only)  is  little 
(Fig.  14(a)),  indicating  that  its  intrinsic 
Iinewidth  is  small.  A  Iinewidth  upper  limit  of 
1 .69  meV  is  determined  for  this  vibrational 
mode.  For  the  CH2  wagging  mode,  a  nonlinear 
least  squares  fit  to  Eq.  (6)  (solid  curve  in  Fig. 
14(b))  gave  an  intrinsic  Iinewidth  of  13.5  ±  2.4 
meV.  The  linewidths  and  their  variation 
throughout  the  molecule  are  potentially  due  to 
inhomogeneous  configuration  of  the  molecular 
constituents,  and  a  more  detailed  understanding 
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may  give  detailed  structural  information  of  these  device  structures. 


Figure.  14.  Line  broadenings  as  function  of  AC  modulation  for  IETS  spectra  obtained  at  4.2  K  for 
(a)  the  Au-S  stretching  mode  and  (b)  the  CH2  wagging  mode.  The  circular  symbols  are 
experimental  FWHMs  and  the  square  symbols  are  theoretical  calculations  considering  modulation 
and  thermal  contributions.  Nonlinear  least  squares  fitting  to  determine  intrinsic  linewidth  is 
shown  as  the  solid  curve  in  (b).  The  intrinsic  linewidths  obtained  for  the  Au-S  stretching  mode  is 
<  1 .69  meV,  and  for  the  CH2  wagging  mode  is  13.5  ±  2.4  meV. 
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Numerous  talk  and  seminars.  No  consultations  or  transitions 


New  discoveries,  patents 


none 


Addendum:  Collaboration  with  Xia,  U.  Washington 


Through  this  DURINT,  we  obtained  samples  from  Xia  et  al,  to  investigate  electrical  and 
optical  effects.  As  a  result  of  this  collaboration,  we  had  a  serendipitous  and  important 
discovery. 

We  have  looking  at  simple  approaches  for  both  coupling  to  propagating  plasmon 
modes  and  their  observation  in  metallic  nanowires.  Specifically,  our  experiment  couples 
to  propagating  modes  in  a  radially  symmetric  nanowire  by  using  one  end  as  a  scattering 
center.  This  scattering  center  creates  an  overlap  between  the  incident  radiation  (focused 
laser)  and  the  propagating  plasmon  mode.  Focused  laser  light  excites  plasmons  that 
propagate  along  the  length  of  silver  nanowires  and  couple  to  free  space  photons,  radiating 
at  the  ends.  Plasmons  are  also  observed  to  couple  between  overlapping  nanowires  and 
fan-out  from  one  wire  into  multiple  nanowires. 

The  wires  used  in  this  study  have  a  mean  diameter  of —100  nm  and  lengths  that 
vary  from  3-20  pm.  The  nanowires  were  synthesized  in  a  mixture  of  ethylene  glycol  (EG) 
and  poly(vinyl  pyrrolidone)  (PVP)[1]  The  wires  were  washed  once  in  acetone  to  remove 
the  EG,  four  times  with  ethanol  to  remove  the  PVP,  and  once  with  water.  The  final 
aqueous  suspension  consisted  of  silver  nanowires,  a  small  fraction  of  silver  nanoparticles, 
and  trace  amounts  of  PVP  and  ethanol.  The  structure  of  the  nanowires  was  characterized 
by  transmission  electron  microscopy  (TEM),  as  seen  in  Figure  1.  The  nanowires  are 
observed  to  have  tapered  ends  (Figure  1  A),  and  are  smooth  to  within  the  resolution  of  the 
TEM  micrographs  (~lnm  ).  Additionally,  a  thin  layer  of  PVP  surfactant  (2-5  nm)  can  be 
resolved  in  many  of  the  micrographs.  An  electron  diffraction  pattern  shows  that  the  bulk 
of  the  wire  is  highly  crystalline  (see  Figure  IB)  and  has  face  centered  cubic  (FCC)  lattice 
symmetry. 


87 


Figure  1.  Structure  of  silver  nanowires.  A)  Transmission  electron  micrograph  (TEM)  of  the  silver  nanowires.  B)  Electron  diffraction 
pattern  of  the  same  nanowire,  showing  FCC  crystalline  structure. 


Aqueous  nanowire  suspensions  were  deposited  on  No.l  Zi  cover  glasses  (Corning 
No.  021 1  zinc  titania  glass)  and  allowed  to  dry  in  open  air.  Dried  nanowires  were  then 
mounted  on  an  inverted  optical  microscope  (Nikon  Eclipse  TE-2000).  A  lOOx  oil 
immersion  objective  (N.A.  =  1 .4)  lens  was  used  to  focus  the  laser  light  and  to  collect  a 
bright  field  image.  The  laser  illumination  was  coupled  into  the  microscope  via  a  dichroic 
mirror  that  selectively  reflects  96%  of  light  at  830nm.  The  microscope  objective  focuses 
the  collimated  laser  light  to  a  diffraction  limited  spot  in  its  focal  plane.  Images  were 
collected  with  either  a  CCD  (Hitachi,  8  bit,  480X640)  or  a  high-speed  CMOS  (Photron 
FastCam  1024-PCI),  which  both  received  roughly  4%  of  the  light  from  the  sample  at  the 
laser  frequency  (see  supporting  information).  To  eliminate  the  possibility  of  evanescent 
waves  propagating  along  the  glass  surface  and  preferentially  scattering  from  the  tips  of 
the  silver  nanowires,  we  immersed  all  samples  with  index-matching  oil  (Nikon  Type  A 
immersion  oil,  n=1.515).  Even  though  scattering  from  the  glass  surface  disappeared,  light 
continued  to  strongly  radiate  from  the  distal  end  of  the  nanowire. 

We  launch  plasmons  by  illuminating  an  end  of  a  single  nanowire  with  a 
diffraction  limited  laser  spot  as  shown  in  Figure  2.  Plasmons  can  be  launched  from  either 
end  of  the  nanowire  (Figures  2A  and  2B);  thus,  plasmon  propagation  is  reversible.  In 
contrast,  plasmon  modes  are  not  observed  to  be  launched  when  the  laser  is  focused  on  the 
mid-section  of  a  smooth  wire  (see  Figure  2C  and  2D).  Since  the  momentum  of  the 
incoming  photon  (kphoton)  is  not  matched  to  that  of  the  propagating  plasmon  (kpiaSmon), 
there  needs  to  be  a  scattering  mechanism  to  provide  an  additional  wavevector  (Akscatter). 

At  the  midsection  of  the  wire,  the  nanowire  is  cylindrically  symmetric  over  the  extent  of 
the  diffraction  limited  spot,  and  therefore  cannot  scatter  in  the  axial  direction.  However, 
this  symmetry  is  broken  at  the  tapered  end  of  the  nanowire  where  light  is  scattered  into 
propagating  axial  plasmon  modes.  A  similar  strategy  has  been  implemented  for  thin  film 
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surface  plasmons  utilizing  gratings  or  dots  [2]  Likewise,  propagating  plasmon  modes 
incident  upon  a  sharp  discontinuity  (e.g.,  the  tapered  end  of  the  nanowire)  can  re-emit  as 
photons. _ 


Figure  2.  Micrographs  showing  the  spatial  sensitivity  of  launching  plasmons.  A) 
Nanowire  with  excitation  at  the  bottom  end.  B)  Same  nanowire  when  excited  from  top 
end.  C)  Nanowire  excited  at  left  end.  D)  Same  nanowire  with  laser  positioned  at  the 
middle  of  the  nanowire.  Notice  that  the  plasmon  is  not  excited  in  this  geometry. 


If  the  symmetry  is  gently  broken  over  longer  length  scales,  plasmon  propagation 
is  unaffected.  This  can  be  seen  in  Figures  2A  and  2B,  where  plasmons  propagate  around 
the  bend  of  the  nanowire  with  no  observed  radiative  loss.  The  smallest  naturally 
occurring  optically  resolvable  bend  found  in  these  nanowires  is  shown  in  Figure  3A. 

This  nanowire,  with  a  radius  of  curvature  of  4pm,  guides  plasmons  with  no  observed 
radiative  loss.  However,  extremely  sharp  bends  in  the  nanowire  will  behave  like  the  end 
of  nanowires,  scattering  propagating  plasmons  into  photons.  This  effect  is  demonstrated 
in  Figure  3B,  where  a  nanowire  shows  emission  at  two  “kinks”  (radius  of  curvature 
below  the  diffraction  limit).  The  kinks  are  spontaneous  defects  that  form  during  the 
growth  process,  and  occur  in  a  small  fraction  of  wires.  Scanning  electron  micrographs 
reveal  that  typical  kinks  are  discontinuities  in  the  wire  direction.  Figure  3C  is 
representative  of  the  structure  of  these  kinks;  i.e.,  sharp  interior  angles  and  flat  exterior 
faces,  with  characteristic  size  ~  lOOnm. 
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Figure  3.  Micrographs  of  plasmon  propagation  in  silver  nanowires  and  emission  (top), 
control  with  no  laser  excitation  (bottom).  Brightest  point  in  image  is  scattering  from 
incident  beam.  A)  A  7pm  wire  (excitation  top)  with  a  radius  of  curvature  of  4pm,  a  wire 
not  radiating  is  in  close  proximity.  Inset  is  a  circle  of  radius  of  4um  for  comparison.  B) 
A  5pm  wire  (excitation  top  right)  with  both  the  opposite  end  and  two  additional  points  of 
high  curvature  that  radiate  (radius  of  curvature  less  than  the  diffraction  limit.)  C)  An 
electron  micrograph  of  a  typical  wire  with  kinks,  with  insets  of  increased  magnification. 


A  nanowire  with  multiple  kinks  can  be  used  to  estimate  the  plasmon  propagation  length. 
Each  kink  is  used  dually  as  a  plasmon  launching  site  and  photon  radiating  site.  We 
sequentially  irradiate  each  kink  and  measure  the  radiated  intensity  from  all  the  other 
kinks,  as  illustrated  in  Figure  4A.  We  find  that  the  radiated  intensity  is  strongly 
correlated  to  the  distance  between  points  along  the  nanowire,  as  shown  in  Figure  4B. 

x 

These  data  are  well  described  by  an  exponential  decay,  I(x )  =  I0e  L ,  with  a  characteristic 
length  of  L=3  ±lpm.  If  one  makes  the  assumptions  that  i)  the  photon-plasmon  coupling 
strengths  are  the  same  at  each  junction,  and  that  ii)  radiative  losses  (at  kinks)  are 
insignificant  compared  to  dissipative  losses  (along  the  nanowire  length),  then  L 
represents  the  plasmon  propagation  length.  Even  when  these  assumptions  are  not  true, 
this  simple  far-field  measurement  provides  a  practical  lower  limit  for  the  propagation 
length. 
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Figure  4.  A)  Micrographs  of  plasmon  propagation  in  a  silver  nanowire  with  multiple 
emission  points  (top)  control  with  no  laser  excitation  (bottom).  Camera  shutter  speed  and 
laser  intensity  were  varied  to  increase  total  dynamic  range.  The  unsaturated  images  of 
the  excitation  and  brightest  emission  points  at  shorter  exposure  times  (high  shutter 
speeds)  are  overlaid  for  clarity.  B)  Semi-log  plot  of  intensity  versus  distance  from 
excitation  source.  The  dashed  line  is  a  fit  of  the  data  to  an  exponential.  Six  distinct 
points  were  probed  (15  unique  combinations).  The  uncertainty  is  arrived  at  by  taking  the 
standard  deviation  of  all  exposures  with  unsaturated  pixels. 


Propagating  modes  are  also  seen  to  couple  between  overlapping  nanowires.  In 
Figure  5  we  present  an  image  and  an  intensity  profile  from  three  overlapping  nanowires, 
which  form  naturally  during  solvent  evaporation.  When  one  of  these  nanowires  is 
excited,  two  observations  are  seen.  First,  visible  radiation  is  emitted  at  the  intersection  of 
the  overlapping  nanowires.  Second,  visible  radiation  is  emitted  at  the  ends  of  nanowires 
that  overlap  the  illuminated  nanowire.  This  is  reminiscent  of  the  interwire  coupling  seen 
in  sub-wavelength  silica  and  semiconductor  nanowire.  While  coupling  for  dielectric 
guides  decreases  with  increasing  diameter  since  less  energy  is  available  in  the  evanescent 
wave,  metallic  nanowire  coupling  should  not  have  this  diameter  dependence  since 
plasmons  remain  confined  to  the  surface. 
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Distance  (um) 


Figure  5.  Group  of  overlaying  nanowires  that  illustrates  inter-wire  plasmon  coupling.  The  excitation  at  the  far  left  nanowire  end 
produces  reemission  at  both  the  nanowire  junctions  and  reemission  from  the  coupled  nanowires  as  well.  The  inset  shows  an  intensity 

line  cut  showing  the  reemission  intensity  profile  along  wire. 


In  conclusion,  we  demonstrate  the  selective  launch  and  propagation  of  plasmons 
along  silver  nanowires  using  a  simple  far-field  excitation  and  detection  method.  We  have 
also  observed  that  these  propagating  plasmon  modes  can  couple  between  adjacent 
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nanowires.  The  phenomena  are  not  specific  to  nanowire  material  or  excitation 
wavelength  -  we  have  observed  similar  results  in  gold  and  copper  nanowires. 
Technological  applications  require  the  positioning  of  nanowires  into  a  deterministic 
network  which  can  be  accomplished  by  a  variety  of  methods,  such  as  flow  alignment17, 
biologically  derived  templates18,  dielectrophoretic  alignment19  or  holographic  optical 
tweezers  20,21  The  current  study  represents  a  new  approach  in  the  observation  and 
manipulation  of  plasmons  in  nanoscale  structures. 

1 .  Sun,  Y.;  Yin,  Y.;  Mayers,  B.  T.;  Herricks,  T.;  Xia,  Y.  Chemistry  of  Materials 
2002, 14,  4736-4745. 

2.  Ditlbacher,  H.;  Krenn,  J.  R.;  Schider,  G.;  Leitner,  A.;  Aussenegg,  F.  R.  Applied 
Physics  Letters  2002,  81,  1762. 
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HOLOGRAPHIC  INTERFERENCE  LITHOGRAPHY 


FOR  PHOTONIC  AND  PHONONIC  CRYSTALS 


Chaitanya  Ullal  and  Edwin  Thomas 


Our  broad  goal  of  this  area  of  the  DURINT  is  the  use  of  holographic  interference  lithography  for 
the  fabrication  of  two-dimensional  periodic  and  three  dimensional  periodic  bicontinuous 
structures.  Previously,  work  in  this  area  had  concentrated  on  the  establishment  of  symmetries,  the 
application  of  these  structures  as  photonic  crystals,  and  as  microlenses.  The  research  work  over 
the  past  review  period  has  built  on  these  efforts  and  also  introduced  some  new  thrust  areas.  The 
recent  work  has  been  in  three  areas: 

1.  Photonic  Crystals: 

a.  Development  and  fabrication  of  architectures  relevant  to  on  chip  photonic 
crystals 

b.  Writing  of  defects  using  2-photon  lithography,  (collaboration  with  SUNY- 
Buffalo) 

c.  Development  of  characterization  for  specific  3d  bicontinuous  structures 
(collaboration  with  WPAFB) 

2.  2D  holographic  patterns  as  phase  masks. 

3.  Phononic  crystals:  development  of  a  complete  tool  set  for  the  study  of  hypersonic 
phononic  crystals 

1.  HOLOGRAPHIC  INTERFERENCE  LITHOGRAPHY  FOR  PHOTONIC 
CRYSTALS 

In  previous  work  in  this  DURINT  we  had  examined  the  effect  of  symmetries  on  the 
photonic  band  gap  for  various  structures.  We  had  found  that  structures  with  simple 
low  Fourier  terms  tend  to  support  gaps.  This  work  subsequently  led  to  a  deeper,  more 
physical  insight,  by  viewing  these  low  Fourier  term  structures  as  possessing 
“dielectric  mirrors”  in  principal  directions.  In  order  to  arrive  at  experimentally 
fabricable  structures  that  support  band  gaps,  we  restricted  ourselves  to  structures  that 
have  sinusoidal  modulations  in  real  space  that  are  along  reciprocal  lattice  vectors  that 
are  the  closest  to  the  origin.  We  were  thus,  interested  in  targeting  three  such 
structures:  viz.  the  level  set  approximations  to  the  Schwartz  Simple  Cubic  (P),  Gyroid 
(G),  and  Diamond  (D).  As  previously  demonstrated  these  structures  are  theoretically 
accessible.  However  from  the  point  of  view  of  their  use  as  structures  to  be  made  on 
chip,  given  the  opacity  of  the  substrate,  we  required  all  the  interfering  beams  to  come 
from  the  same  half  space.  In  order  to  negotiate  this  consideration  we  derived 
approximations  to  these  structures  that  retained  some  of  the  sinusoidal  band  gap 
supporting  terms  while  also  being  amenable  to  fabrication.  We  demonstrated  that  the 
level  set  approximation  to  the  P  structure,  the  3-FCC-term  and  3-BCC-term  structures 
were  three  such  structures.  Importantly  the  3-FCC-term  structure  displayed  two  gaps 
and  the  P  structure  is  scalable  in  terms  of  unit  cell  size,  for  any  fixed  writing 
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wavelength.  These  two  structures  were  then  successfully  fabricated  in  an  epoxy 
based  negative  photoresist. 

A  more  recent  area  of  thrust  in  terms  of  photonic  crystals  has  been  the  controlled 
introduction  of  defects  via  two-photon  lithography.  In  order  to  introduce  defects  into 
the  photonic  crystal  structures  being  written  by  interference  lithography  we  add  a  step 
to  the  fabrication  procedure.  Once  the  photoresist  has  been  exposed  to  the 
interference  pattern  the  sample  is  taken  over  to  a  two-photon  direct  writing  setup. 
Defects  are  written  into  the  photoresist  by  scanning  the  focal  point  of  an  infra  red 
laser  beam  into  the  volume  of  the  photoresist.  The  photoresist  is  then  subjected  to  the 
previously  established  processing  steps,  viz.  post  exposure  baking,  developing  and 
supercritical  drying.  The  net  pattern  is  the  sum  of  both  the  defects  and  the  large  area 
single  crystalline  interference  lithography  pattern.  Preliminary,  results  are  shown  in 
Figure  1.  This  work  is  being  done  as  a  collaborative  effort  between  the  Thomas 
group  at  MIT  and  the  Prasad  Group  at  SUNY-Buffalo. _ _ _ 
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An  important  area  of  interest  in  the  fabrication  of  three  dimensional  structures  via 
holographic  interference  lithography  is  the  development  of  a  non  destructive 
characterization  technique  in  the  estimation  of  volume  fraction.  This  is  particularly 
important  in  the  case  of  photonic  crystals  since  the  photonic  band  gaps  are  strongly 
dependent  on  the  volume  fraction  of  the  structures.  We  are  currently  developing  a 
technique  based  on  diffraction  to  determine  volume  fraction  of  particular  structure. 
This  involves  the  transfer  of  interference  patterns  into  h-pdlc  materials  and  the 
subsequent  characterization  of  these  materials  via  diffraction.  This  work  is  being  done 
in  collaboration  with  the  Bunning  group  at  WPAFB. 

2.  2D  holographic  lithographic  microlens  patterns  as  phase  masks. 

In  previous  work  under  this  DURINT  we  had  examined  the  potential  for  using  2d 
holographic  patterns  to  create  biomimetic  microlenses.  These  lenses  were  further  shown 
to  behave  as  microfluidic  elements  with  integrated  pores.  We  also  extended  the  use  of 
these  biomimetic  patterns  by  demonstrating  their  application  as  multipattem  photomasks; 
that  is,  by  using  the  same  photomask  and  simply  adjusting  (i)  the  illumination  dose,  (ii) 
the  distance  between  the  mask  and  the  photoresist  film,  and  (iii)  the  tone  of  photoresist, 
we  were  able  to  create  a  variety  of  different  microscale  patterns  with  controlled  sizes, 
geometries,  and  symmetries  that  originated  from  the  lenses,  clear  windows,  or  their 
combination  (See  Figure  2).  Light  field  calculations  were  made  and  found  to  be  close  to 
the  experimental  results. 
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3.  Hypersonic  Phononic  crystals. 

A  new  thrust  area  with  this  DURINT  in  the  area  of  holographic  structures  has  been 
the  study  of  hypersonic  phononic  crystals.  To  date  the  focus  of  the  research  in  this 
area  has  been  in  the  application  of  the  holographic  structures  as  photonic  crystals. 
This  is  based  on  the  fact  that  we  have  periodicities  that  are  on  the  length  scale  of 
visible  light.  We  proposed  and  demonstrated  the  use  of  hypersonic  phononic  crystals 
in  the  control  of  emission  and  propagation  of  high  frequency  phonons.  We  reported 
the  fabrication  of  high  quality,  single  crystalline  hypersonic  crystals  using 
interference  lithography  and  showed  that  direct  measurement  of  their  phononic  band 
structure  (See  Figure  3)  is  possible  with  Brillouin  light  scattering.  Numerical 
calculations  were  employed  to  explain  the  nature  of  the  observed  propagation  modes. 
This  work  lays  the  foundation  for  experimental  studies  of  hypersonic  crystals  and, 
more  generally,  phonon-dependent  processes  in  nanostructures. 


triangles,  glass  mode;  open  triangles,  Bragg  mode;  solid  circles,  phononic  crystal  modes;  solid  lines, 
theoretical  quasilongitudinal  modes;  dotted  lines,  theoretical  mixed  modes. 
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POLYMERIC  PHOTONIC  CRYSTALS 
FOR  LOW-THRESHOLD  PHOTONIC  BANDGAP  LASING 


Jongseung  Yoon,  Wonmok  Lee  and  Edwin  L.  Thomas 

The  objective  of  this  area  of  the  DURINT  is  the  use  of  polymeric  photonic  bandgap 
structures  for  fabricating  low-threshold  photonic  bandgap  laser  devices.  We  pursued  the 
development  of  defect-mode  one-dimensional  (ID)  photonic  bandgap  lasers  by  two 
different  approaches:  (1)  block  copolymer  self-assembly  and  (2)  polymer/inorganic 
nanoparticle  hybrid  layer-by-layer  assembly. 


1.  Optically  Pumped  Surface-Emitting  Lasing  using  Self-Assembled  Block 
Copolymer  Distributed  Bragg  Reflectors 

With  an  appropriate  domain  size  that  is  comparable  to  optical  wavelengths,  block 
copolymers  can  be  used  to  create  periodic  dielectric  structures  which  can  interact  with 
visible  light,  i.e.  photonic  crystals.  Over  the  past  years,  block  copolymers  have  been 
considered  as  a  unique  materials  platform  for  fabricating  large-area,  well-ordered 
photonic  bandgap  structures.  One  potential  application  of  such  self-assembled  polymeric 
photonic  crystals  is  to  use  them  as  a  resonator  in  a  photonic  microcavity  to  provide 
spectrally-selective  feedback  for  lasing.  In  this  work,  we  fabricated  a  thin  film  organic 
laser  cavity  using  a  block  copolymer  based  one-dimensional  ID  photonic  crystal. 
Polymeric  distributed  Bragg  reflectors  (DBRs)  were  prepared  through  the  self-assembly 
of  a  lamellar-forming  poly(styrene-b-isoprene)  (PS-b-PI)  diblock  copolymer  having  a  ID 
photonic  stop  band  overlapping  with  the  fluorescence  spectrum  of  a  gain  medium. 
Optically  pumped  surface-emitting  lasing  was  obtained  using  polymethylmethacrylate 
(PMMA)  doped  with  l,4-di-(2-methylstyryl)benzene  (Bis-MSB)  as  an  organic  gain 
medium  and  the  polymeric  self-assembled  DBR  as  a  spectral-band  selective  feedback 
elements  (See  Figures  4,  5,  6).  This  block  copolymer  based  photonic  structure  opens  the 
possibility  for  creating  all-organic,  flexible,  and  self-assembled  active  photonic  devices 
with  fast  and  low-cost  processing,  and  spectral  tunability  by  various  external  fields. 


98 


BCP  DBR 


Gain  medium 
Cover  glass 
Slide  glass 


Figure  4.  Schematic  of  the  block  copolymer  based  laser  cavity,  comprised  of  a  gain  medium,  Bis-MSB  and 
PMMA,  enclosed  between  two  block  copolymer  based  distributed  Bragg  reflectors.  The  thicknesses  of  PS 
and  PI  domains  are  estimated  values  at  10  wt  %  solvent  concentration  based  on  the  peak  position  and 
FWHM  of  the  reflectivity  spectrum. 
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2.  Defect-mode  Mirrorless  Lasing  in  Dye-doped  Organic/Inorganic  Hybrid  ID 
Photonic  Crystal 


We  developed  a  novel  dye-doped  organic/inorganic  hybrid  ID  photonic  crystal 
containing  a  dye-doped  defect  layer  for  defect-mode  photonic  bandgap  lasing.  The 
multilayer  laser  structure  consists  of  alternating  layers  of  titania  (TiCL)  nanoparticles  and 
polymethylmethacrylate  (PMMA)  with  an  active  emission  layer  of  organic  dyes  (DCM) 
in  PMMA  (See  Figures  7,8,9).  Low  threshold  lasing  has  been  demonstrated  at  a  single 
defect-mode  wavelength  of  the  ID  photonic  bandgap  structure  resulting  from  the 
inhibited  density  of  states  of  photons  within  the  stop  band  and  the  enhanced  rates  of 
spontaneous  emission  at  the  localized  resonant  defect  mode.  This  work  has  been  done  as 
a  collaborative  effort  between  the  Thomas  group  at  MIT  and  the  Prasad  Group  at  SUNY- 
Buffalo. 
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Figure  7.  Schematic  of  dye-doped  defect-mode  ID  photonic  crystal,  glass-(PMMA-Ti02)15- 
(DCM/PMM  A)-(Ti02-PMM  A) '  5-air. 
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Figure  8.  The  lasing  spectrum  obtained  at  a  pump  power  of  1  mW,  above  the  lasing  threshold.  The  two 
small  peaks  beside  the  lasing  line  at  582  nm  correspond  to  the  excitation  light  (532  nm)  and  to  the  low 
frequency  defect  mode  (at  620  nm,  below  its  lasing  threshold  for  this  pump  power). 
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Figure  9.  Emission  intensity  and  line-width  (the  full  width  at  half  maximum,  FWHM)  at  the  lasing 
wavelength  (582  nm)  as  a  function  of  pump  power.  The  behavior  of  the  intensity  and  the  FWHM  clearly 
demonstrate  a  threshold  for  lasing  around  0.6  mW  pump  power  (12  pJ  pulse  energy). 
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